


Examples of waves breaking from a storm in the North Atlantic in December 1993 in which the winds were gusting up to ap- 
proximately 50-60 knots and wave heights of up to 12-15m were reported. a: Large scale breaking waves. b: Smaller-scale 
waves breaking on longer waves, c: Short strongly wind-forced breaking waves. (Photographs by E. Terrill and W. K. Melville, 
Scripps Institute of Oceanography) 
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This Naval Research Reviews on the topic of nonlinear ocean waves is pub- 
lished in two companion issues, Two/1996 (pages 1-28) and Three/1996 
(pages 29-72). The introduction which puts the six articles into context 

and provides an historical perspective appears at the beginning of Two/1996. 
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2D Complex Wavelets for 
Analysis of Unsteady, Wind- 
Generated Surface Waves 


G.R.Spedding 


Departments of Aerospace & Mechanical Engineering 


University of Southern California Los Angeles, CA 


Abstract 


Wavelets and their transforms (WT) have numerous and 
diverse applications, and almost as many shapes and sizes. 
This article concentrates solely on their utility as a quantitative 
measure of spatially-localised spectral amplitudes for 2D data 
analysis. The physical problem concerns the time-dependent 
growth of wind- driven surface waves, starting from a flat 
surface. Although the problem lacks the full complexity of real 
ocean conditions, it is argued that because it is precisely 
defined, a combined experimental and theoretical program can 
advance our physical understanding, that may itself be pro- 
jected onto field data. Certainly, the advent of nonlinearity in 
the waveficld is marked by intermittency and 3D structure, and 
the WT is an appropriate basis for analysis. The wind-wave 
problem will be used as an example where the simultaneous, 
quantitative measurement of amplitude and phase modulations 
in an unsteady wave field leads to estimates of quantities such 
as the local dispersion relation that cannot easily be measured 
otherwise. Furthermore, the WT-derived results can be com- 
bined with simple analytical models to measure quantities that 
are both of great interest to the oceanographer, and that have 
traditionally been difficult to estimate. Finally, suggestions are 
made for future work and take a surprising direction. 
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What is the problem? 


Why study small-scale wind waves? 


For the purposes of this article, small-scale waves are 
those having wavelengths in the centimeter range. Both gravi- 
tational and surface tension effects are significant and so the 
physical problem is an interesting balance between the two. 
From a practical point of view, there are two reasons for 
interest. Capillary-gravity (cg) waves are the wave-like distur- 
bances that first grow when a calm water surface is disturbed 
by wind. Later in time, the energy- containing scales are much 
larger, but at any instant, it is through smaller-scale motions 
such as cg waves that the instantaneous surface response 
mediates the fluxes of momentum and energy between the 
atmosphere and ocean. The second practical motivation is that 
remote sensing devices such as SAR are sensitive to scattering 
of electromagnetic radiation primarily in the shorter wave- 
lengths characterised by cg waves. It is only by modulation of 
the shorter wavelengths that larger scale features are observ- 
able. Generally it is sufficient to model this effect as a ‘rough- 
ness’ , but statistical roughness models and predictions will not 
behave realistically unless they are based on sound physical 
principles. As we shall see later, the cg wavefield can be 





extremely complex in structure, and further understanding 
leads inevitably to studies that fall squarely into the ‘Nonlinear 
Dynamics’ theme of this issue. 


Why study the unsteady problem? 


It is a mild understatement to observe that fluctuations of 
the local wind stress (in both space and time) occur over the 
ocean surface. While steady-state response models may work 
satisfactorily in some statistical sense, the transient response 
and physical understanding of how one state changes to an- 
other, together with the possible interactions between a time- 
varying component with another stationary one, are still of 
great interest. Furthermore, while the linearised models predict 
the initial stages of the surface response quite well, they cannot 
predict the subsequent evolution of low wavenumber (long 
wavelength) modes, 3D pattern instabilities, and wave-wave 
interactions that follow. In other words, there still is a problem, 
and it has not yet been solved. 

There is a second, more abstract, line of reasoning that 
also supports such a line of enquiry. One of the principal 
achievements of nonlinear dynamical systems theory has been 
to demonstrate how complex phenomena can be generated by 
simple, deterministic systems, with few degrees of freedom. 
Hopes have been raised that tasks of apparently daunting 
complexity may yet yield to analytical, predictive methods. In 
a similar vein, starting from the very evident fact that the real 
sea surface state is extremely complicated, one searches for the 
simplest possible model that can nevertheless exhibit suffi- 
ciently rich dynamics that it can account for the most important 
physical mechanisms. 

Starting from scratch, we require a system that: (i) is 
clearly defined; (ii) can be modeled analytically, or numeri- 
cally; (iii) can be tested experimentally, and (iv) can be gener- 
alised to situations of practical interest. Condition (iv) is rarely 
satisfied unless (ii) and (iii) involve a degree of understanding 
based on physical principles, and the failure of empirical 
matching of model parameters and coefficients when applied 
outside of the domain for which they have been carefully tuned 
is a rather familiar problem in oceanography. 

The spectrum of analytical models for the air-water inter- 
face is bounded on the one hand by small amplitude, two-di- 
mensional solutions of the Rayleigh equation for steady 
mean flows with boundaries at x = +e. Correspondingly, 
the Blasius solution to the incompressible, flat- plate bound- 
ary-layer equations for a specified origin in x gives the 
fetch-limited result. Between the two limiting cases, we are 
interested in the time-dependent problem, in the absence of 
external boundaries. 

One might question the relevance of laminar model flows 
when the ocean applications are almost always turbulent, in 
both fluids. However, in the initial transient for a wind profile 
starting from rest, even if the flow itself is not perfectly 
laminar, the laminar instability modes could still be predomi- 


nant in selecting the most energetic wave modes. In this 
respect, the laboratory experiment, incorporating some degree 
of turbulence in the growing profiles, provides a bridge be- 
tween model and reality. 


Why bother with laboratory experiments? 


In the previous article, Norden Huang has considered the 
general applicability of wave-tank experiments to real ocean 
conditions. It is worth reiterating the basic philosophy behind 
our laboratory work and its relation to the overall problem 
outlined in the previous paragraphs. In contrast to field data, it 
is only in the relatively controlled conditions of the laboratory 
that parameters can be varied at will, independently, and re- 
peatedly. It is not necessarily significant that conditions do not 
exactly match those of the open ocean, because the experi- 
ments can still provide non-trivial tests of analytical or con- 
ceptual models that are themselves of use over a broader range 
of parameters. The support or refutation of such models can 
still provide information about their ability to predict events in 
the real world. If and when the model fails, it can be extended 
or thrown away; in either event, the conditions under which it 
fails are now known 

In the remainder of this article, a closely-linked theoreti- 
cal and experimental research program will be described, 
where wavelet analysis will be used as a tool for identifying 
current limitations and pointing to future requirements. 


The contribution of wavelet- 
based methods 


Wavelets as atomic units 


At the time of writing, there are 11 texts on wavelets in 
our university bookstore, and some are quite good. For a 
concise and readable introduction, it is hard to beat the Algo- 
rithms and Applications book by Meyer'. A more comprehen- 
sive and mathematical overview is presented in the 
introductory volume of the series on wavelets edited by Chui’, 
and geophysical applications of varying quality can be found 
in a later volume’ in the same series. Farge‘ has written a fine 
review on wavelets and turbulence, and further collections of 
theoretical advances and applications are accumulating.>” 

One of the great strengths of the field is in the mathemati- 
cal depth and rigour of much of the derivation and analysis, 
which has revealed connections with hitherto unrelated tech- 
niques, from Litthewood-Paley analysis of the 1930’s to mod- 
ern quadrature mirror signal processing algorithms and spline 
functional analysis. Nevertheless, a zero-th order degree of 
sophistication in understanding wavelet-based signal analysis 
will suffice for our purposes. On this level, wavelets are simple 
oscillatory functions that have compact support, and thus can 
be used for local space-scale decompositions. The original 
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function, f(x) (or f(2)), is expanded as a function of an addi- 
tional scale variable, a 


fix) Ff (ax). (1) 


The decomposition is effected by a convolution of f(x) 
with multiple copies of a single wavelet function, w (x), that 
are systematically shifted and dilated to cover a range of 
characteristic scales, a, and positions, b 

Provided certain admissibility criteria are met w (x) must 
have finite energy, zero mean...consult refs.(1,2,4)), then the 
choice of w (x) is somewhat arbitrary. It is possible to make 
the strong assertion that the physical world is much more likely 
to be composed of collections of wavelet-type, compactly 
supported functions, rather than, say, infinitely- extended os- 
cillatory functions as in Fourier series. It follows that it is only 
reasonable to decompose real world observations (and com- 
putations) onto such atomic units. Farge’, for example, draws 
attention to the fact that a decomposition of a complex fluid 
flow into its elementary coherent structures (if such an object 
exists), is better approximated by a functional base with com- 
pact support than the usual Fourier base. Alternatively, one 
might search for basis functions whose properties are extracted 
from the data, or from known solutions to the particular 
governing equations. There will likely be future progress in 
this direction, but it is in fact unnecessary to adhere so strictly 
to the atomic construct. 


Wavelets as tools for quantitative 
description 


The keyword, of course, is ‘construct’, and one is free to 
construct a wavelet, or any other, basis to represent the world 
in many different ways. For example, one of the principal 
advantages of Fourier representations is the considerable 
mathematical convenience in computation of PDE’s, having 
little or nothing to do with any implied physical basis. Simi- 
larly, provided one is aware of the properties of the decompo- 
sition and its role in determining (even distorting) the 
interpretation of the world, one can still use a wavelet basis for 
descriptive purposes. This is not all the same thing as stating 
that wavelet-atoms actually exist, or that they have any special 
dynamical significance. Hence, we simply state the following: 


1. The objective is to make quantitative measurements and 
interpretations of local energy densities in a 2D, time-varying 
surface wave field. 

2. The most straightforward way to do this is through 
continuous, complex-valued wavelets. 


The definition of the wavelet transform (WT), described 
in words above, can be statd more precisely: 


Fab) =| fix).y fz * ae (2) 


The transform is a sum of the convolution of the translates 
and dilates of w*(x) with f(x) (the * denotes the complex 





Figure 1. 


Is the world composed of wavepackets like this? The Morlet 
wavelet is a smooth, continuous wave function, modulated by 
a Gaussian envelope. Real and imaginary parts are shown in 
continuous and dotted lines. Whether or not it is a true atomic 
unit of the physical world, it can still be a convenient function 


basis for general description. 
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Figure 2. 


The 2D wavelet function Arc is defined most conveniently in 
Fourier space where for a given scale, a, it can be seen to 
interrogate an arc that covers all wavevectors of magnitude 
Ikl, regardless of their orientation. 
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conjugate). In one dimension, a useful choice of w(x) is the 
Morlet® wavelet. 


w(x) = eo" - er), (3) 


As Figure 1 demonstrates, it is acomplex-valued function, 
and the 1/2 phase shift between real and imaginary parts allows 
the modulus of its coefficients to be free from oscillations due 
to changes in phase in an otherwise constant amplitude signal. 
The wavelet phase indicates the phase in the instantaneous 
frequency, or wavenumber, associated with a particular scale, 
a. 

k, is constant (typically) so there is a fixed number of 
cycles in an analysing wavelet, regardless of a. Consequently, 
the spatial resolution improves with increasing wavenumber, 
or decreasing a. The way in which it does so varies automat- 
ically with a, and the balance between resolution in wavenum- 
ber space and in physical space is in some sense an optimum 
for the Morlet wavelet, within the overall constraints of Heis- 
enberg’s uncertainty principle. 

The WT conserves energy, and a form of the Parseval 
relation can be written, 


fimmPdr=2-f firanipani 22, @ 


2 , 
a 





Figure 3. 


The time evolution of alongwind surface slope following a 
change in mean wind, U1o, from Om/s to 6.8m/s. The change 
in wind speed is approximately linear over 5 seconds. The time 
series shown here begins at t = 5.3s, i.e. after Uio has reached 
its steady value. The field of view in each image is approxi- 
mately 14 x 11.5 cm. 





where Cy is a normalisation constant that depends on the form 
of y (x ). The right hand side is an integral of a function that 
is proportional to the spatial distribution of energy density per 
unit scale. It is possible to integrate contributions from all x to 
arrive at a global power spectrum, 





Figure 4. 


Experimentally-measured initial exponential growth rate, B, vs. 
wavenumber, k,(symbols connected by solid lines), for U1o 
from 2.8 - 10.2 m/s. The solid symbols are from experiments 
of Kawai, and the dashed lines are results from Kawai'’s Orr- 
Sommerfeld model integration for parameters that closely 
matched those experiments (Uo is given in brackets). At 
moderate windspeeds, the agreement is good, but is very 
sensitive to assumptions about how the boundary layer grows, 
and so remains somewhat empirical, because there is no 
analytical model for the boundary layer growth. A simple 
diffusion model does not work at all, and a pessimist might well 
argue that the agreement simply signifies that both parties 
have done similar experiments. As selected, the best match 
occurs at moderate windspeeds, but the extremes are not so 
well covered. 
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yI Ss. (5) 


E(a) = c,! flab 


E(a) is a Gaussian-smoothed version of a power spectrum that 
could equally well have been obtained from Fourier trans- 
forms. Of course, in this case we are not obliged to perform 
the summation over x, and the point is to note that the energy 
densities are precise, quantitative measures of localised spec- 
tral components (provided that there are no artifacts from the 
mutual interaction of overlapping, non-orthogonal wavelet 
coefficients). 

The extension to two dimensions of this basic approach 
is fairly straightforward and details can be found in ref. 9. 
As before, we require a complex wavelet function, and for 
real-valued signals (2D surface slope measurements, for 
example) it is sufficient to cover the half-plane in wavenum- 
ber space. Now that yw=w(x,y), the wavelet can be designed 
so that it can be either directional-specific, or non-direc- 
tional. A non-directional wavelet responds equally to 
wavevectors at all angular orientations, and in two dimen- 
sions such a function covers an arc that sweeps through 180° 
in (k,, k,} (Fig. 2). 

Arc and its close relatives form the core of the data 
analysis, giving a measure of instantaneous, wavenumber 
energy densities and wave phase as a function of {x,y}. 


An integrated analytical/ 
experimental program 


Theory 


The time-varying solution for perturbations to the in- 
terface between two viscous fluids with specified shear- 
flow profiles can be solved numerically by integrating 
coupled Orr- Sommerfeld equations for air and water. Both 
Valenzuela’® and Kawai'' have done that, two decades or so 
ago. Exponential growth rates for cg waves were predicted 
(and measured, in the case of Kawai) and compared reason- 
ably well with microwave backscatter measurements of 
Larson & Wright.'* This body of work represents, still, the 
state-of-the-art in the numerical analysis of the linearised 
problem. 

The stability of piecewise linear shear profiles in air and 
water was recently examined by Caponi et al.'* One can take 
their steady-state model and integrate forwards in time by 
allowing a boundary layer thickness to grow, thus generating 
a quasi-steady solution. The sensitivity of the solution to the 
assumed growth of the boundary layer thickness will tell us 
much about the behaviour of the numerical solutions, and the 
directions for further improvements. 
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Figure 5. 


The modulus of the wavelet transform, |WTI,,, at the initially 
most amplified scale (k,=4.0 rad/cm, 4,=1.6cm). The colour 
bar mapping is constant over the series, which covers the 
exponential growth, saturation, and decline in energy densities 
around this scale. Although the wavenumber and growth rates 
are quite well predicted, as described in Figure 4, the distribution 
is far from uniform {x,y}. 








Experiment 


Experimentally, we investigate the growth of disturbances 
on an initially calm water surface when the wind starts from 
rest, and at sufficient fetch so that upstream boundary condi- 
tions can be ignored (at least in the water). The configuration 
is a reasonably close experimental approximation of the theo- 
retical approaches discussed above. 

The measurements are made with a modified form of the 
Imaging Slope Gauge (originally described by Jahne & Rie- 
mer'*), which uses surface refraction of intensity-graduated 
incident light to derive the instantaneous along- or cross-wind 
slope over the field of view. s,(x,y,t) or s,(x,y,t) can be recon- 
structed according to the direction of the light intensity gradi- 
ent. The surface slope, s,,, is directly related to the wave 
energy, and the WT time series, 5, ,(x,y,a,t) contains informa- 
tion about the local energy densities, and about the local 
wavenumber- frequency relation. Model predictions can thus 
be carefully checked. For the first time, the unsteady, local 
dispersion relation can be measured, and related to the 2D, 
instantaneous structure of the wavefield. 


Matching the two 


Neither existing models, nor practicable experiments can 
provide complete and generalisable information about the 
time-varying, 3D, wave and velocity fields. It will emerge that 
the theoretical models are very sensitive to, and experiments 





mostly ignorant of, the detailed shape and time-evolution of 
the water and air boundary layers. However, in combination, 
mutual consistency checks and judicious matching of empiri- 
cal data with analytical relations can lead to new methods for 
measuring and/or inferring the values of physical variables 
(surface stress, boundary layer height), not only from labora- 
tory data, but also from remote sensing measurements. Two 
examples of this kind of reasoning will be given. In both cases, 
the glue that holds the pieces together is the Arc wavelet 
transform. 


The time evolution of capillary- 
gravity waves 


The surface wave growth for a mean wind, U,, = 6.8 m/s 
(measured at 10 cm above the calm water surface) is shown in 
Figure 3. The time series proceeds from left to right, top to 
bottom, and is for wind forcing from the right. The dimension- 
less time Tt = U,of/X, where X is the fetch, ranges from 3 to 5.5, 
so the upstream origin is always felt in the air flow, even at the 
early times. The group velocity of the initial water waves, 
however, is much smaller. At 50 cm/s, for example, 24 seconds 
are required to travel the fetch distance of 12m. 

It is clear that, even at early times, the observed wavefield 
is never perfectly two-dimensional. The grey-scale encoding 





Figure 6. 


IWT |a-1 for same time series as Figure 5. Maximum growth of 
the lower wavenumber mode (k-1 = 2.2 rad/em, h-1 = 2.9 cm) 
occurs around the middle of the third row, just when the peak 
of the initially most amplified mode in Figure 5 begins to 
saturate and decay. The boxes marked ‘k;’ are located at the 
energy density peaks and correlate well with peaks at the 
higher wavenumber mode (k1=7.9 rad/cm, 49=0.8 cm; not 
shown). 








of the slope magnitude gives the initial disturbances an organ- 
pipe appearance, but they always have Y-shaped ‘defects’, or 
‘dislocations’, where successive waves become entangled. In 
fact, it was the presence of the defect patterns that first drew 
our attention to the wind-wave problem, as we wondered 
whether there may be similarities with other defect structures 
in fluid turbulence.!*! 

Even though the waves are not perfectly two-dimensional, 
calculation of the global power spectral densities and their rate 
of change with time shows that the most amplified wavenum- 
bers, and their initial exponential growth rates, can be quite 
well predicted by linear theory (Figure 4). That is true for the 
first five images (the top row) in Figure 3. During the course 
of about 1.25s, the initially most-amplified disturbances have 
reached their maximum amplitude, and in the top right panel, 
the waves look as two-dimensional as they ever will be. Shortly 
afterwards, the surface patterns become much more compli- 
cated as the initial modes saturate and further growth occurs 
in sidebands of both higher and lower wavenumber. By the end 
of the series longer wavelength disturbances predominate. A 
convincing physical explanation for this well-known phe- 
nomenon has never been given. Note, however, that it does 
appear to be a duration-limited event, since information propa- 
gated by water waves from upstream boundaries will not yet 
have reached the measurement window. Further insights can 
be derived from wavelet analysis. 


Wavelet-based analysis 


Local energy density distributions 


The remarks concerning the global evolution of instability 
modes are based on Fourier analysis. The 2D Fourier transform 
can provide a perfectly good description of the average prop- 
erties of the surface slope, and has very good spectral resolu- 
tion. Having identified the wavenumbers of interest, it is now 
possible to use the localised properties of the wavelet trans- 
form to relate the spectral densities to the ‘structure of the 
wavefield. 

The spatial distribution of growth of the initially most 
amplified mode at U,, = 6.8 m/s can be seen in Figure 5. The 
energy density distributions are quite patchy, and remain so. 
By analogy with defects in mixing layers, where the defects 
are the sites for the first (and subsequent) subharmonic transi- 
tions, it is reasonable to ask whether the globally-measurable 
sideband instabilities are located at, or near, defect sites in the 
surface wave field. The open circles denote the location of 
defects, identified by eye in the equivalent time series of 
surface slope. Again, by analogy with the mixing layer results, 
the local energy density at the fundamental scale, ay, would be 
a minimum here, while the sideband modes would peak at the 
same location. At first glance, the circles do seem to be 
consistentlv beside, rather than on, the highest energy densities 
in ay. However, Figure 6, which is the same plot for the lower 
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Figure 7. 


The waveprint. The sign of the phase of the WT at ao allows 
defect sites to be located with little ambiguity, enabling auto- 
mated identification and tracking algorithms to work reliably. 
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wavenumber sideband, a,, shows that the circles are not 
associated with local maxima at this scale. It can also be found 
that the other sideband, characterised by wavelet scales a,, is 
correlated quite well with maxima in @_,, and so, also not with 
a». The high and low wavenumber modes correlate well with 
each other, but only weakly with minima in the initial mode. 
Initially, the defects do not appear to play a significant role in 
the development of sidebands. Although one cannot expect 
spatial correlations to be maintained in a dispersive wavefield, 
the variation in phase speed with wavenumber is rather weak 
at this particular range of scales, and the wave pattern defects 
or dislocations do not appear to be involved in the growth of 
sidebands. 








Pattern defects and the waveprint 


The correlation of pattern defects with various postulated 
transitions and instabilities could be made more precise if a 
reliable method for identifying the defect sites were available, 
and we recently found one. 

Recall that it is specifically defects at a particular scale 
that we would like to investigate, and the WT can be employed 
as a kind of filter so that only contributions from this scale are 
counted. Images such as Figures 5 & 6 are cleariy of no use, 
but the phase of the WT is. Figure 7 is a binarised image of the 
phase of the WT at scale a,. Not only can one confirm the 
defect locations that were originally identified by eye (circles), 
but it is also possible to follow their time evolution. None of 
the defect topologies in the circles remains fixed over time; 
instead connections with preceding or following waves are 
initiated and broken. In a similar spirit to the description of 
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fluid flows by critical points,'* the wavefield topology can be 
characterised by the type and location of its defects, together 
with the distribution of the highest amplitude wavelet coeffi- 
cients. Just as fingerprints can be categorised by their defects, 
so the wavefield can be described by its waveprint. One can 
look for particular defect structures in an appropriate analytical 
model, and Paul Newton at USC has made a start in describing 
a new family of defect solutions to the nonlinear Schrodinger 
equation.'* 





Figure 8. 


How to calculate the instantaneous, local dispersion relation 
w/(k, x, y) from a time series of complex wavelet transforms. 
Given the surface slope distribution, Sx(x,y) (top row), at two 
closely-spaced times, t1, tz, the wavelet transform of each can 
be separated into modulus and phase components (rows 2 & 
3). For all scales having significant energy densities, the flow 
field that maps the isophase lines from their position at t1 to 
their position at te can be computed from an optical flow 
technique. Note that frequencies higher than the apparent 
Nyquist limit can be found because the structure in the phase 
field allows large displacements to be computed in a global 
optimisation solution. Finally, only measurements that occur 
where IWTI > |WTlihresh are retained. 
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u(x,y) for IWTI > thresh 








Combined analytical/ 
experimental models for 
prediction 


Estimation of local dispersion relation 


The WT gives a measure of the spatially-localised energy 
densities as a function of wavenumber, and since the measure- 
ments are also well-resolved in time, the local frequencies can 
be calculated from @ = 0 / dt. This is possible because the 
complex WT also gives the instantaneous phase distribution, 
(x,y), associated with each scale (as shown in the waveprint). 
Figure 8 shows how the calculation proceeds. The flowfield, 
u(x,y), that maps $(x,y,t,) to (x,y,,) at all energetic scales, a 
(associated with physical wavenumber, k), is used to compute 
the local frequency as. 


@ (x,y) =| k (xy)li7tx,y). (6) 


There is no comparable time-dependent, 2D theory that 
predicts such a detailed dispersion relation, but values over all 
{x,y} can be collected together and compared. Figure 9 is an 
example, where the solid curves have been drawn over an 
arbitrary range of k from the result of Caponi et al'* for the 
intrinsic @(k) of small perturbations on a surface between two 
piecewise linear shear profiles. The experimental results differ 
from the theoretical curves, which is not altogether surprising 
since the wavelengths (and frequencies) that actually have the 
largest amplification rate in practice will depend on the time 
history of the growth in boundary layer thickness, D(t), and 
surface drift velocity, U,(t). In fact, from the difference be- 
tween the measured @(k) and the solid curve, either the surface 
drift or the boundary layer thickness can be inferred. This is 
an example of a combined experimental/theoretical modeling 
for simultaneous measurement of a variable whose estimation 
could be of great importance in the field. 


Estimation of surface stress from lab, 
or field data 


Now let us suppose that U, is known (from existing 
experimental data or field measurements. The maximum slope 
thickness, D, can be written as the ratio of the mean wind, U, 
to the gradient of dU/dy at the water surface 


The friction velocity, u, , is a measure of the wind stress 
exerted at the water surface, which also varies with dU/dy: 


(8) 


where 1 and p are the air viscosity and density, respectively. 
Combining the two expressions, the friction velocity can be 
related to D, 


(9) 


where v=1/p is the kinematic viscosity. An estimate of D is 
therefore equivalent to an estimate of the friction velocity. To 
see how this relationship can be incorporated into the overall 





Figure 9. 


An example comparison of measured @(k) with theoretical 
predictions in solid curves. Each symbol marks the average 
value of the collection of points at a particular wavenumber 
where IWTI exceeds a threshold value. Each spatial location 
contributes a dot, and the collection of dots forms a vertical 
line whose length shows the variation in measured @ for a given 
k. The first two timesteps cover the regime that is adequately 
described by linear theory. Aftert=4.95, the scatter increases 
and by t=5.23, there are no wavelet coefficients that exceed 
the threshold value at the originally most-amplified wavenum- 
ber. 
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framework, see Figure 10. At the heart of the scheme is the 
combined theoretical model for w(k) together with experimen- 
tal measurements of the same quantity. Given a measured 
dispersion relation, and a selection criterion for choosing one 
particular curve, such as the maximum amplification rate, then 
the difference between measured and predicted @(k) curves 
can be attributed to the surface drift current, U,,. The corrected 
ratio of the water and air velocities are returned to the theoreti- 
cal model, and upon convergence, a predicted value of D=D, 
can be used to measure u,. It is noteworthy that the measure- 
ments do not have to come from laboratory data. In principle, 
a time series of SAR images could provide the same informa- 
tion, because the WT can give estimates of w(k) from these 
data also. 

The two examples above demonstrate how the newly-ac- 
quired ability to measure instantaneous amplitude and phase 
information for local spectral components can lead to novel 
and potentially powerful measurement techniques that have 
application in both laboratory and ocean data. 


Concluding remarks 


The wind-wave startup problem has been explored as an 
example of the application of wavelet-based measurement and 
analysis techniques to unsteady, inhomogeneous data. It is the 
extension to these domains that brings us closer to realistic 
ocean conditions, and clearly into nonlinear regimes that have 
just begun to be investigated. The advanced measurement 
techniques have application to both laboratory and field con- 
ditions. Sie!) 

With regard to the particular wind-wave problem, it is 
obvious that the wavefield is always quite complex in struc- 
ture, but the dynamical significance of the complexity is not 
equally obvious. In attempting to distinguish between struc- 
ture that is due to intrinsic instability of the water surface 
deformation, or that merely reflects imprecision in specifica- 
tion of the initial forcing, we need a smaller-scale, more easily 
controllable experiment. The form of the mean wind profile, 
and its time history, must be specified. Then, an experimental 
program with modulations in both along- and cross-wind 
directions to generate particular instability modes would be a 
wind-forced equivalent of existing systematic investigations 
of liquid surface instabilities.2°” 

Paradoxically, the conclusion that future experiments 
ought to be in smaller, more precisely- controlled containers, 
moves still further from true-blue oceanography, where 
Reynolds numbers are large, and many simultaneous proc- 
esses compete and/or interact with each other. Recalling the 
introduction however, it is likely that the basic instabilities and 
physical mechanisms will only be isolated and identified with 
a more selective, and precisely specified experiment. We hope 
to report on the resulis of similar experiments in the future. 
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Figure 10. 


A flow diagram for the procedure for estimating wind friction 
velocity from measured inputs, @, k and Ug. The inputs may 
come from either laboratory or field data. The subscripts w and 
a refer to quantities in water and air, respectively, and sub- 
scripts m and p identify measured and predicted values. 
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Introduction 


Wind-generated surface gravity waves (figure 1, frequen- 
cies between approximately 0.05 and 0.5 Hz) are the energy 
source for nearshore fluid motions having a range of spatial 
and temporal scales. In shallow water (depths less than about 
10 m) near-resonant nonlinear interactions result in significant 
cross-spectral energy transfers and strong evolution of wave 
shapes over distances of only a few wavelengths. As waves 
shoal, energy is nonlinearly transferred from the wind waves 
(about 0.1 Hz in figure 2) that dominate the wave field outside 
the surf zone (736 cm depth in figure 2) to both higher and 
lower frequencies. Sum interactions excite phase coupled mo- 
tions at relatively high frequencies, resulting in the sharp crests 
and steep forward faces characteristic of near-breaking waves 
and bores (153 cm depth in figure 2). Difference interactions 
are believed to generate lower frequency (nominally 0.005- 
0.05 Hz) infragravity waves (“surf beat”) which dominate 
current and sea-surface elevation fluctuations in the swash 
zone where wave breaking has substantially reduced the wind- 
wave energy (0 cm depth in figure 2). 

These nonlinearly excited motions may play a central role 
in nearshore sediment transport. Phase coupling between wind 
waves and higher frequency harmonics causes asymmetrical 
near-bottom horizontal velocities and accelerations widely 
believed important to cross-shore sediment transport (eg, 


Bowen 1980, Bailard 1981, Bailard & Inman 1981, Hanes & 
Huntley 1986, Doering & Bowen 1988, Elgar et al. 1988). 
Strong infragravity motions are the dominant sediment sus- 
pension mechanism in the inner surf zone during storms (eg, 
Holman et al. 1978, Wright et al. 1982, Beach & Sternberg 
1988). 

The nonlinear evolution of ocean surface gravity waves 
propagating through the nearshore to the beach observed at 
several field sites will be reviewed here. First, sum interactions 
that transfer energy to higher frequency motions will be con- 
sidered, followed by a discussion of the excitation of lower 
frequency motions by difference interactions. 


Nonlinear Sum Interactions 


Well seaward of the beach, in depths comparable to or 
greater than the wavelength (\klh > 1, where k is a typical 
wavenumber and / the water depth), ocean waves are charac- 
terized by strong frequency dispersion and broad directional 
distributions. Changes in the wave field caused by gradual 
variations in the water depth are well described by linear 
propagation models (eg, Longuet-Higgins 1957, Collins 1972, 
Berkhoff 1972, Radder 1979, Le Me’haute’ & Wang 1982, 
Kirby 1986a, 1986b, Kirby & Dalrymple 1986, Panchang et 
al. 1990). At second order in weakly nonlinear theory (Phillips 
1960, Hasselmann 1962), forced secondary motions arise 
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Figure 7. 


Two breaking waves (in about 2 and 5 m depths) at the U.S. 
Army Corps of Engineers’ Field Research Facility pier near 
Duck, NC. The deck is about 7.5 m above mean sea level. The 
surf zone extends more than 1 km offshore. (From Elgar et al. 
1995.) 








which can be interpreted as small, but significant corrections 
to the underlying linear wave field. For example, interactions 
between directionally opposing wind waves excite long wave- 
length secondary waves that dominate high frequency (eg, 
greater than about twice the spectral peak frequency) pressure 
fluctuations at the seafloor (figure 3), which in tum drive 
seafloor microseisms (Longuet-Higgins 1950). Energy levels 
and wavenumbers of high frequency, secondary pressure fluc- 
tuations observed at the seafloor in relatively deep water (IklA 
=3) are consistent with the predictions of weakly nonlinear 
theory (Herbers & Guza 1991, 1992, 1994). Resonances be- 
tween quartets of waves occur at the next higher order, result- 
ing in slow cross-spectral energy transfers. Although nonlinear 
energy exchanges are small on wavelength scales, the fre- 
quency-directional spectrum can be substantially modified 
over hundreds of wavelengths (Phillips 1960, Hasselmann 
1962, and others). As the depth decreases (Iklh < 1), secondary 
wave effects are enhanced and skewed wave shapes develop 
with sharp crests and broad, flat troughs. The nonlinearity of 
wave-induced pressures and orbital velocities is readily de- 
tected in intermediate water depths and accurately predicted 
by second-order nonlinear theory for a wide range of wave 
conditions (Herbers et al. 1992). 

In very shallow water (Iklh < 1, typically near or within 
the surf zone),'waves are approximately nondispersive, direc- 
tional distributions are narrow (owing to refraction during 
shoaling), and strong nonlinearities drive relatively rapid spec- 
tral evolution. Models based on the nondispersive, nonlinear 
shallow water equations for unidirectional waves (Hibberd & 
Peregrine 1979) predict many important features of bore 
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propagation and runup (Kobayashi et al. 1989, Raubenheimer 
et al. 1995, and references therein). 

The transition region between dispersive deep and non- 
dispersive shallow water regimes is the shoaling region, char- 
acterized by weakening dispersion, narrowing directional 
spread, and substantial nonlinearity. Waves propagating 
through the shoaling region evolve significantly in several 
(rather than hundreds of) wavelengths, with nonlinear interac- 
tions driving phase-locked harmonics that cause the pitched- 
forward shape characteristic of nearly breaking and broken 
waves (Freilich & Guza 1984, Elgar & Guza 1985a, 1985b, 
Elgar et al. 1990, Kaihatu & Kirby 1995, and many others). In 
particular, nonlinear interactions occur among triads of waves 
with frequencies f, and wavenumbers k; obeying the interac- 
tion rules (Armstrong et al. 1962, Bretherton 1964): 


k, +k, —k,=0 (1b) 


For sum interactions, the primary wave components 1 and 2 
each obey the lowest-order (linear) dispersion relation 





Figure 1A. 


Elevation versus time (arbitrary units) of a sawtooth wave 
shape (dashed line, skewness = 0, asymmetry = 2.3) and its 
Hilbert transform (solid line, skewness = 2.3 and asymmetry 
= 0.). The time series have identical power and bicoherence 
spectra, but different biphases. (From Elgar et al. 1990.) 
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Figure 2. 


Near-bottom pressure time series (top) and power spectra 
(bottom) well outside the surf zone (736 cm depth), in the 
breaking region (153 cm depth), and on the beach face 
(swash excursions, 0 cm depth). 
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and force a new component 3 with the sum frequency and the 
vector-sum wavenumber. In the nonresonant case where com- 
ponent 3 does not obey (1c), there is a small energy transfer to 
secondary “bound waves” (Phillips 1960, Hasselmann 1962). 
On the other hand, in the resonant case where the sum compo- 
nent exactly satisfies the dispersion relation (1c), continuous 
energy transfers occur between the three modes and compo- 
nent 3 grows at the expense of components 1 and 2. Armstrong 
et al. (1962) showed that significant energy also can be ex- 
changed between the three waves of the triad if the sum 
component nearly satisfies the dispersion relation (“near-reso- 
nant” interactions). 


The evolution of a wave field dominated by nearly nor- 
mally incident, narrow-band swell (peak frequency f = 0.06 
Hz, significant wave height H, = 60 cm in 9 m depth) observed 
on a monotonically sloping (slope ~ 0.05) beach near Santa 
Barbara, CA is shown in figures 4-6. In 9 m water depth (figure 
4a) there is a small spectral peak at the first harmonic (f= 0.12 
Hz) of the swell. As the waves shoal, the power in the first and 
higher harmonic (f ~ 0.12, 0.18, 0.24, 0.30 Hz) peaks increases 
(figures 4b-d, upper panels). The growth of these harmonics is 
not predicted by any linear theory, but is accurately modeled 
(Freilich & Guza 1984) by near-resonant sum interactions 
based on the nonlinear Boussinesq equations (Peregrine 1967). 





Figure 3. 


Bottom-pressure spectral density versus frequency at 4-min 
intervals (times and line types are given in the legend) imme- 
diately following a sudden shift in wind direction. Note the large 
increases in secondary-pressure spectral levels at double- 
wind wave frequencies (about 0.35-0.70 Hz). The pressure 
gage was located on the seafloor in 13 m depth near Duck, 
NC. Pressure spectral have been converted to equivalent 
spectra of vertical water column displacements. (From Her- 
bers & Guza 1994.) 
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Figure 4. 


Power spectra (crr*/Hz) and contours of bicoherence for 
narrow band swell. The power spectra (bars indicate 95% 
confidence ievels) are immediately above the corresponding 
bicoherence spectra. Contours indicate quadratic coupling 
among the triad of waves with frequencies f1, fo, and f; + fa. 
The minimum bicoherence contour level is b = 0.1, with 
additional contours every 0.05. There are 310 degrees of 
freedom and the 95% significance level for zero bicoherence 
is b = 0.14. Depths are (a) h = 9.0, (b) 3.9, (c) 2.7, and (d) 1.3 
m. The data were observed on Leadbetter Beach, Santa 
Barbara, CA. (From Elgar & Guza 1985b.) 
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Below each power spectrum in figure 4 is the correspond- 
ing bicoherence spectrum, b ( f, f, ) which detects phase cou- 
pling between triads of waves with frequencies 
ff, and f, +f, (see Appendix for details on bispectral analy- 
sis). If the 3 waves in the triad are independent of each other, 
as would occur in a Gaussian (linear) sea, b* ( f, f, )=0. On 
the other hand, b” (f, f, ) = 1.0 if all the energy with frequency 
f, +f, is phase coupled to the waves with frequencies 
f, and f,. In a natural wave field interactions occur across a 
broad range of frequencies so that 0<- (f,,f,)<1, and 
b (f,,f) can be interpreted as a relative measure of nonlinear 
coupling between waves with frequencies f, ,f,, and f, +f,. 

At the seaward edge of the shoaling region (water depth 
h = 9 m, figure 4a) there is weak coupling between the swell 
(f, =f, = 0.06 Hz) and its harmonic (f,= 0.12 Hz). As the waves 
shoal, the excitation of phase-coupled harmonics is clearly 
visible in the bicoherence spectra. In 4 m depth (figure 4b), the 
bicoherence indicates enhanced coupling between the swell 
and its first two harmonics (triads with frequencies (0.06, 0.06, 
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Figure 5. 


Biphase versus depth for selected frequency pairs (with 
nonzero bicoherence) for the narrow band swell. Asterisk, (f, 
f); circle, (f, 2f); plus, (f, 3f); triangle, (2f, 2f); diamond, (2f, 3f), 
where f = 0.06 Hz corresponds to the power-spectral-peak 
frequency. (From Elgar & Guza 1985b.) 














0.12) and (0.06, 0.12, 0.18) Hz). In 1.3 m depth (figure 4d) the 
unbroken swell is coupled to its first through fifth harmonics 
(contours with f= 0.06 Hz), and the harmonic with frequency 
0.12 Hz is coupled to even higher harmonics (contours with f, 
= 0.12 Hz). Although these bispectral calculations indicate 
which frequencies are nonlinearly coupled, and not the direc- 
tion of energy flow (ie, which frequencies are receiving en- 
ergy), the sequence of power spectra (figure 4) shows that 
energy is predominantly transferred to higher frequencies. 

Along with increased bicoherence during shoaling (figure 
4), there is substantial biphase evolution (figure 5). The 
biphase is a statistical measure of the shapes of the waves (see 
Appendix). For example, the biphase of a Stokes wave (a wave 
with tall, narrow crests and broad, flat toughs) is zero, while 
the biphase of a sawtouth is — 2/2. The bicoherence and 
biphase evolution (figures 4 and 5) show that Stokes-like 
waves in 9 m depth (the biphase of the most energetic triad 
f, =f, = 0.06 Hz is ~ 0) evolved to sawtoothed-shaped waves 
with steep front faces and gently sloping rear faces (biphases 
of all triads are ~ — 2/2) in 1 m depth. These observed shape 
changes are also evident in the time series that show wave 
profiles evolving from slightly peaked shapes in 9 m depth to 
the pitched-forward shapes of nearly breaking waves in shal- 
lower water (figure 6). 

The observed evolution of a locally generated broadband 
wave field (figure 7, also with H= 60 cm in 9 m depth) is 
markedly different. The power spectra and bispectra are fea- 
tureless, showing no evidence of the prominent harmonic 
peaks observed with the narrow-band swell (compare figures 
4 and 7). The bicoherences evolved from near-zero values in 
9 m depth to low (but statistically significant) values spread 
over the entire wind-wave frequency range in 0.9 m depth 








Sea-surface elevation versus time for a short section (about 
90 s) of the narrow bank swell data set. (a) h = 8.7, (b) 2.9, (c) 
2.4, and (d) 1.0 m. (From Elgar & Guza 1985b.) 
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(figure 7). Strong nonlinearity occurs as in the narrow swell 
case, but it is masked in the power spectra because the inter- 
actions are distributed over a wide range of frequencies. 

Despite the radically different evolution of power and 
bicoherence spectra for the narrow and broadband wave fields, 
the evolution of biphase and wave shape is similar. In both 
cases, the biphases evolved from near zero in 9 m depth to 
— 1/2 just before the waves broke, and the waves had steep 
front faces and flat rear slopes in very shallow water (compare 
figures 6 and 8). 

A bulk statistical characterization of an average wave 
shape is provided by third moments of the sea-surface eleva- 
tion time serie’. Skewness, the normalized third moment of the 
time series, indicates the degree of peakedness of the waves 
(eg, the lack of symmetry about a horizontal plane caused by 
differences between the shapes of wave crests and troughs). 
Asymmetry, another normalized third moment of the time 
series, indicates the degree of pitched-forwardness of the 
waves (eg, the lack of symmetry about a vertical plane caused 


by differences in the shapes of leading and trailing wave faces). 
Skewness and asymmetry are equal to normalized sums of the 
real and imaginary parts of the bispectrum, respectively. Thus, 
bispectra (see Appendix) indicate which frequency triads con- 
tribute to total third moments (skewness and asymmetry) of 
measured time series, similar to the way in which the standard 
frequency spectrum shows which frequencies contribute to the 
total second moment (variance). Although the bicoherence 
spectra for the two wave fields are quite different, the evolution 
of the bulk wave skewness and asymmetry is similar. Skewness 
is relatively low in 9 m depth (figure 9), increases to a maxi- 
mum in about 2 m depth, and then decreases to near zero in 
shallower water. Asymmetry is near zero in 9 m depth, but 
increases as the waves shoal. 

The observations shown in figures 1-9 were obtained on 
monotonically and moderately sloping beaches where wave 
breaking limits the spatial extent over which nondissipative, 
nonlinear evolution can occur. Significantly stronger nonlinear 
evolution was observed over a nearly horizontal (flat), 2 m 
deep beach section extending about 80 m from the toe of the 
foreshore seaward to a small sand bar (figure 10) during the 
recent ONR funded Duck94 nearshore field experiment. In 
particular, nearly breaking (“feathering” of the wave crests was 
observed) narrow-band swell (H, = 80 cm) that was strongly 
pitched forward near the sand bar (asymmetry ~ 0.5 near 





Figure 7. 


Power spectra and contours of bicoherence for a broadband 
wave field, in a format identical to figure 3. Depths are (a) h = 8.7, 
(b) 3.5, (c) 22, and (d) 0.9 m. The increase in high frequency 
energy (f > 0.2 Hz) is an artifact of not correcting the deepest 
measurement (made with a bottom-mounted pressure gage) for 
depth attenuation. (From Elgar & Guza 1985b.) 
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cross-shore coordinate x = 280 m, figure 10) became unpitched 
during further propagation in constant depth. At the shoreward 
end of the flat section (x = 180 m, figure 10) the asymmetry 
was nearly zero (and even slightly negative), in contrast to 
moderately sloping beaches where strong positive asymmetry 
(pitched forwardness) is a precursor to breaking and waves 
remain pitched forward across the surf zone. The changes in 
the wave spectrum during propagation across the flat section 
of the beach were dramatic (figure 11). In 8 m depth (cross- 
shore coordinate x = 884 m in figure 10) the dominant wave 
period was about 11 s (figure 11). In 2 m depth, at the seaward 
edge of the flat section (x = 240 m), small crests were begin- 
ning to emerge in the troughs of the 11 s swell. These reso- 
nantly excited waves continue to grow as the nonbreaking 
waves propagate over the flat beach section, and by its shore- 
ward edge there are twice as many wave crests as offshore 
(compare figures 1la and llc). Although the dramatic fre- 
quency doubling shown in figure 11 was rare because narrow- 
band swell occurred infrequently, nonbreaking waves with a 
range of spectral shapes had similar evolution of skewness and 
asymmetry over the flat beach section. 





Figure 8. 


Sea-surface elevation versus time for a short section (about 
90 s) of the broadband wave field. (a) h = 8.4, (b) 2.5, (c) 1.9, 
and (d) 0.9 m. (From Elgar & Guza 1985b.) 
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Figure 9. 


Skewness and asymmetry versus depth. Asterisk, narrow 
band swell; circle, broadband wave field; solid curves are 
skewness and dashed curves are asymmetry. The data were 
band-pass filtered between 0.04 and 0.4 Hz. (From Elgar & 
Guza 1985b.) 
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There is considerable interest in waves in the infragravity 
frequency band (0.005 < f < 0.05 Hz) on the beach and 
continental shelf. In addition to sometimes dominating wave 
runup at the shoreline (figure 2), infragravity waves can reso- 
nantly excite ocean structures and small harbors. Nonlinear 
difference frequency interactions between a pair of swell com- 
ponents with slightly different frequencies f and f + 8f drive an 
infragravity frequency (Sf) secondary (forced) wave. The mag- 
nitude of this secondary wave, which is phase-locked to groups 
of swell, depends on the directional properties of the swell 
waves (Hasselmann 1962). Predictions of forced wave energy 
levels, based on weakly nonlinear theory and the frequency- 
directional spectrum of swell (estimated with a 24-element 
array of pressure sensors in 13 m depth near Duck, NC) are 
consistently lower than observed infragravity levels (figure 








Figure 10. 


(a) Cross-shore transect of instruments deployed near Duck, 
NC. The symbols represent colocated sonar altimeters (to 
determine the bottom location), bottom-mounted pressure 
sensors, and bidirectional current meters. (b) Skewness and 
asymmetry versus cross-shore location for the narrow band 
swell observed to undergo frequency doubling. Asterisks, 
skewness; squares, asymmetry. The data were band pass 
filtered between 0.05 and 0.3 Hz. 
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waves are predominantly radiated from shore and refractively 
trapped on the shelf. Infragravity wave generation appears to 
be sensitive to the local beach bathymetry, whereas the refrac- 
tive trapping depends on the regional shelf topography. The net 
result of this sensitivity is that free infragravity energy levels 
vary widely between sites in the same water depth with similar 
swell conditions (Herbers et al. 1995b). 

Longuet-Higgins & Stewart (1962) suggested that while 
incident swell are dissipated by breaking in very shallow water, 
the associated forced secondary infragravity waves are some- 
how released as free waves and reflect from the beach. Con- 
sistent with this mechanism, infragravity motions in the surf 
zone have been shown to be phase coupled to groups of wind 
waves seaward of the breaking region (Guza et al. 1984, Elgar 
& Guza 1985b, List 1992). Additional support for the hypothe- 
sis of Longuet-Higgins and Stewart (1962) is provided by the 
observed strong dependence of the ratio of upcoast to down- 
coast directed infragravity wave energy fluxes on the same 
ratio in the swell frequency band (figure 15). Note however 
that the ratio of infragravity energy fluxes is usually much 
closer to 1 (ie, comparable upcoast and downcoast fluxes) than 
is the ratio of swell energy fluxes, suggesting that upcoast (or 
downcoast) propagating swell drives both upcoast and down- 
coast propagating infragravity waves. Herbers et al. (1995a) 
show that secondary infragravity waves (excited by difference 
interactions of swell) exhibit this unusual directional behavior, 





12a). The observed infragravity motions are a mixture of 
forced waves phase locked to the local groups of swell and 
freely propagating (uncoupled) waves. Using bispectral analy- 
sis to separate free from forced infragravity waves, Herbers et 
al. (1994) demonstrate that second-order nonlinear theory 
accurately predicts the observed forced infragravity energy 
levels (figure 12b). However, forced waves usually contribute 
only a small fraction of the total infragravity energy (0.1 to 
30% at the 13 m depth Duck, NC site). Analysis of data from 
many shelf sites (figures 13 and 14) shows that free waves are 
usually the dominant infragravity motion on the shelf (Herbers 
et al. 1995b). However, the relative contribution of forced 
wave energy to the infragravity band is highly variable and 
generally increases with both increasing swell energy and 
decreasing depth (figure 14). During high energy conditions 
(H, > 4 m) forced wave energy levels may exceed free wave 
levels even in 200 m depth (figure 14c). 

Strong correlations between energy levels of infragravity 
waves and swell observed at all sites (figure 13) indicate that 
free infragravity waves are driven by swell, as suggested by 
Munk (1949) and Tucker (1950) nearly 50 years ago. However, 
the details of the infragravity wave generation process and the 
subsequent propagation of infragravity waves across the shelf 
are still poorly understood. Observed directional properties 
and cross-shore energy variations show that free infragravity 





Figure 11. 


Sea-surface elevation versus time for the narrow-band, fre- 
quency-doubling swell. (a) cross-shore location x = 884, (b) x 
= 240 (the seaward edge of the flat section), and (c) x = 160 
m (the shoreward edge of the flat section). There are twice as 
many wave crests at x = 160 than at x = 240 m. (d) Energy 
density versus frequency at the same locations. 
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Figure 12. 


(a) Predicted forced versus observed total (free plus forced) 
infragravity energy and (b) predicted versus observed forced 


infragravity energy in 13 m depth near Duck, NC. (From 
Herbers et al. 1994.) 
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and that the observed up to downcoast infragravity energy flux 
ratio agrees well (figure 16) with predictions based on the 
hypothesis that forced secondary waves are released during the 
breaking process. Although further work is needed to assess 
the importance of wave breaking to infragravity wave genera- 
tion and damping, these results lend strong support to the basic 
hypothesis that infragravity waves are driven by nonlinear 
difference interactions between pairs of swell components. 


Summary 


Near-resonant nonlinear triad interactions result in sig- 
nificant cross-spectral energy transfers and phase modifica- 
tions to ocean surface gravity waves over distances of a few 
wavelengths in water depths less than about 10 m. As waves 
shoal, energy is nonlinearly transferred from wind waves 
(frequencies between about 0.05 and 0.2 Hz) that dominate the 
wave field outside the surf zone to both lower and higher 
frequencies. Nonlinear sum-frequency interactions among 
wind waves excite phase coupled motions at higher frequen- 
cies, resulting in the sharp crests and steep forward faces of 
near-breaking waves and bores. Nonlinear difference-fre- 
quency interactions generate lower frequency (nominally 
0.005-0.05 Hz) infragravity waves which sometimes dominate 
current and sea-surface elevation fluctuations in the runup 


where wave breaking has substantially reduced the wind-wave 
energy. 


Appendix: Bispectral Analysis 


Since their introduction more than thirty years ago 
(Hasselmann et al. 1963), bispectral techniques, which isolate 
nonlinear interactions between the Fourier components of a 
time series, have been used to study many systems (see Nikias 
& Raghuveer 1987 and Elgar & Chandran 1993 for reviews). 
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Figure 13. 


(top) Cross-shelf depth variations and pressure gage locations 
(asterisks) at several field sites (Chesapeake Light is off the 
Virginia coast; Duck is off the North Carolina Coast; Ventura, 
Redondo, and Santa Rosa are in the Southern California Bight; 
Harvest Platform is off the California coast near Pt. Conception; 
Oahu is off the Hawaii coast.) (center) Average spectral shapes 
(the average of spectra normalized by the variance in the range 
0.0005-0.14 Hz). (bottom) Observed correlation between loga- 
rithms of spectral levels in narrow (0.0005 Hz wide) frequency 
bands with total swell energy (0.05 - 0.14 Hz) as a function of 
frequency. (From Herbers et al. 1995b.) 
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The discrete bispectrum, appropriate for discretely sampled 
data, is (Haubrich 1965, Kim & Powers 1979) 


Bf fy) =E [Ay Ay Ay > f,) ’ (Al) 

where the A,, are complex Fourier coefficients at frequency 
f,, the subscript n is a frequency (modal) index, asterisk 
indicates complex conjugation, and E [-] is the expected-value, 





or average, operator. Similarly, the power spectruin is defined 
here as 


P(f)=ZEIAyAy). 


(A2) 


It is convenient to recast the bispectrum into its normalized 
magnitude and phase, called the squared bicoherence and 
biphase, given respectively by (Kim & Powers 1979) 


IBA fa)P 


b (f fr) = ] ry? 
EMA; A, PETA; 4, PI 


Figure 14. 





(A3) 





(a) Forced infragravity energy, (b) free infragravity energy, and 
(c) the ratio of forced to free infragravity energy versus swell 
energy. The data are from 8 m depth at Duck, NC (triangles, 
upper clouds), 30 m depth at Ventura, CA (asterisks, middle 
clouds), and 204 m depth at Harvest Platform, CA (squares, 
lower clouds). Least-squares-fit curves to the logarithms of the 
observed energies are denoted by solid lines. Dashed lines 
labeled 1 and 2 indicate a linear and quadratic dependence, 
respectively. (From Herbers et al. 1995b.) 
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Figure 15. 


The ratio AS, Fiwn of the upcoast to the downcoast compo- 
nent of the alongshore energy flux in the infragravity frequency 


band versus the analogous flux ratio Fig" FERS), in the swell 
frequency band. The dashed line corresponds to equal ratios. 
The data were observed in 13 m depth near Duck, NC. (From 
Herbers et al. 1995a.) 
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For a three-wave system, Kim & Powers (1979) show that 
b? (ff) represents the fraction of power at frequency f, + f, 
owing to quadratic coupling of the three modes 
(f, f,,and f, +f,). The biphase is related to the shape (in a 
statistical sense) of the time series (Masuda & Kuo 1981, Elgar 
& Guza 1985b, Elgar 1987). 

Consider the time series 





(A 4) 


N()=cos(f,t + 8,) + cos(fat + 62) + Cos(fyt +83) (AS) 


where f, =f, +f). If the 3 Fourier phases are independent and 
randomly distributed in [0,2 1], then the time series will have 
Gaussian statistics and b(f, f,)=0. If 8,and 6, are randomly 
distributed in [0,2 2] and 6, = 8, + 6, + a, where © is a con- 
stant, then although each of the 3 Fourier phases is randomly 


distributed, the phase relationship between the 3 Fourier 
components is not random. The _  biphase 
BF, fa) = 8; — (0, + 8.) =a and bff, f,)=1.0. Time series 
given by (A5) have identical power spectra regardless of the 
presence or absence of phase coupling, but their bispectra 
differ. Phase coupling occurs in systems where quadratic in- 
teractions between the components at frequencies f, andf, 
excite motions at their sum frequency f,. 
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Figure 16. 


Predicted versus observed ratios of upcoast to downcoast 
alongshore energy fluxes in the infragravity frequency band 
AS, Figwn for 20 cases selected to span a wide range of ratios. 
(From Herbers et al. 1995a.) 
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Time series with biphase B(f,f,)=0 have nonzero 
skewness (the mean of the cube of the time series, normalized 
by the variance), and are characterized by wave profiles that 
have sharp peaks, and broad, flat troughs (an asymmetry with 
respect to a horizontal axis, figure Al). The contribution of 
individual frequency triads to the total third moment is given 
by the bispectrum, similar to the way in which the contribution 
of individual frequency bands to the variance is given by the 
power spectrum (Hasselmann et al. 1963). The skewness is the 
sum of the real parts of the bispectrum normalized by the 
variance. For a time series with B(f, f,) = —1/2, the real part of 
the bispectrum and the skewness are zero, but the imaginary 
part of the bispectrum is nonzero and the time series is asym- 
metrical with respect to a vertical axis (eg, dashed line in figure 
Al), and is characterized by fore-aft asymmetric wave forms 
(eg, the steep front faces and gently sloping rear faces of nearly 
breaking ocean waves). Asymmetry, a quantity analogous to 
skewness, is the sum of the imaginary parts of the bispectrum 
normalized by the variance ( Elgar & Guza 1985b, Elgar 1987). 

In general, the Fourier components of a nonlinear random 
process are partiaily coupled, as opposed to the complete 
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coupling discussed above. If the Fourier component of the sum 
frequency (f,) of the example time series in (A5) is given by 


cos(fzt + 63) = Beos(fyt +03.) + V1 -F cos(f,t + 63,) (A6) 


where 6,. = 6, + 6, + a and6, is randomly distributed, then at 
frequency f, the time series contains both a coupled component 
(with phase 6,.) and an uncoupled component (with phase 
@,,). The bicoherence squared b7(f, f,) =? is equal to the 
fraction of the energy at f, that is phase coupled to f, f,. The 
power spectrum is identical to that for the other cases discussed 
above. 

Many observed time series have large background noise 
levels, but random noise (and nonphase-coupled signals) do not 
contribute to the bispectrum. However, the detection of weak 
nonlinear phase coupling in the presence of significant noise is 
limited by the statistical uncertainty in bispectral estimates ob- 
tained from finite length data records. For example, even a truly 
Gaussian process will yield a nonzero bispectrum estimate if a 
finite-length time series is examined. Beginning with Brillinger 
(1965), Rosenblatt & Van Ness (1965), and Brillinger & Rosen- 
blatt (1967a, b) there have been many studies of the statistics of 
estimates of bispectra. These are not reviewed here. Haubrich 
(1965), Hinich & Clay (1968), Kim & Powers (1979), Elgar & 
Sebert (1989), Chandran & Elgar (1991), and others discuss the 
statistics of estimates of bicoherence and biphase. References to 
extensions to higher-order (eg, trispectra) spectra can be found in 
Elgar & Chandran (1993). 
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Solitons and the Inverse 
Scattering Transform: 
Nonlinear Fourier 


Analysis for Shallow 
Water Wave Trains 


A. R. Osborne 
Istituto di Fisica Generale dell’ Universita 
Via Pietro Giuria 1, Torino 10125, Italy 


Fourier analysis is the mathematical technique which 
allows one to spectrally represent a measured time series as a 
linear superposition of sine waves. The sine waves may be 
thought of as a set of linear basis functions. However, the 
fundamental Eulerian equations of motion for surface water 
wave dynamics are fully nonlinear and sine waves are therefore 
not the ideal set of basis functions. Is there a way to Fourier 
analyze wave data in order to be completely consistent with 
the nonlinear Euler equation dynamics? While this problem 
has yet to be completely resolved, herein I discuss how, for 
shallow water wave motion, linear Fourier analysis may be 
nonlinearly generalized using the inverse scattering transform 
(IST): In this formulation the basis functions are the well- 
known classical cnoidal waves of Korteweg and deVries’. I 
discuss specifically how a nonlinear wave train may be repre- 
sented in terms of a linear superposition of cnoidal waves plus 
their mutual nonlinear interactions. This remarkable simplifi- 
cation of the IST formulation has been somewhat of a surprise, 


particularly because of the complicated nature of periodic 
inverse scattering theory”©. The strength of Osborne’ was that 
such a simplification is possible. The strength of the present 
paper is that such a simplification has already been made a 
practical reality. I synthesize several nonlinear wave trains to 
illustrate the method and then nonlinearly Fourier analyze 
shallow water wave data from the Adriatic Sea. 


Introduction 


Important developments in the understanding of nonlinear 
wave propagation have occurred during the past 30 years. Of 
particular relevance to the present work is the discovery of a 
class of physical systems which are called integrable non- 
linear, partial differential wave equations*''. Equations of this 
type are said to be exactly solvable in terms of a relatively new 
method of mathematical physics known as the inverse scatter- 
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ing transform (1ST). A number of these integrable equations 
are applicable to oceanographic situations. These include: 

(1) The Korteweg-deVries (KdV) equation, which de- 
scribes unidirectional shallow water wave propagation. 

(2) The Kadomtsev-Petviashvili (KP) equation, which 
describes two-dimensional shallow water wave propagation. 

(3) The nonlinear Schroedinger (NLS) equation, which 
describes the one-dimensional dynamics of the wave envelope 
function in both shallow and deep water. 


A major focus of the present paper is to give a non-mathe- 
matical introduction to certain aspects of inverse scattering 
theory and to address shallow water wave dynamics from the 
point of view of the KdV equation. Of particular importance 
is the discovery, about 20 years ago, of the theory of periodic 
and almost periodic solutions of integrable equations for 
which the mathematical formulation is called the periodic 
inverse scattering transform*''. This work revealed remark- 
able relationships among the Riemann theory of Abclian func- 
tions, the spectral theory of Schroedinger operators and 
algebraic geometry. Exploitation of periodic IST has yielded 
rather general periodic solutions of large classes of nonlinear 
wave equations. The technique may be cast in terms of a kind 
of nonlinear Fourier analysis. It is this latter property which 
is of prime importance for the study of nonlinear, shallow 
water wave data. 

In spite of rather marvellous successes with regard to the 
theory, it goes without saying that daunting mathematical 
difficulties have damped practical applications of the results 
for over 20 years. The study of periodic inverse scattering 
theory is not normally part of the graduate curricula in most 
applied mathematics and physics departments, much less in 
the applied sciences. It therefore should not be surprising that 
applications of IST to natural wave motions have not been 
ubiquitous. The present paper is an attempt to provide a 
simple, physical perspective of some of the major results 
leading to time series analysis of shallow water waves”: !??, 
Not withstanding the level of mathematical difficulty, which 
is not presented here, the physical results turn out to be 
remarkably simple. 


Sine Waves and Linear 
Fourier Analysis 


Linear Fourier analysis allows an arbitrary wave train to 
be represented as a linear superposition of sie waves of 
particular amplitudes, frequencies and phases. The sine wave 
is a travelling wave solution of a linear wave equation and the 
Fourier transform is a spectral solution of the equation (see 
Fig. 1). An example is the following linear shallow water 
equation with dispersion: 


N,+con,+Bn,,,=0 (1) 
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Here 1 (x,t ) is the free surface elevation as a function of 
space x and time f. The coefficients c»=Vgh B=c,h*/6, 


are constants depending upon the depth / and gravitational 
acceleration g. The simplest periodic solution to (1) is a sine 
wave: 


T (x,t) = 9 sin ( kyx — Wot + Po ) (2) 





Figure 1. 


An example of linear Fourier analysis. In (a) is the Fourier 
spectrum, i.e. the amplitudes of the sine waves in the spectrum 
as a function of wave number. In (b) are the sine waves 
themselves; their linear superposition gives the synthesized 
wave train shown at the bottom of the panel. 
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Figure 2. 


Examples of cnoidal waves. In (a) is a solitary wave or soliton, 
in (b) a Stokes wave and in (c) is a sine wave. 
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Here 1) » is the wave amplitude, X, is the wave number, Wp is 


the frequency and @, is the phase. The Fourier approach 
allows one to generally solve (1) for periodic boundary con- 


ditions, for all space and time, by a linear superposition of 
the solutions (2): 


N 
Sn (xt) = > n, sin ( k,x- Ot + ®,) . & 


n=1 


The sine waves in the Fourier spectrum all have unique 
amplitudes, wave numbers, frequencies and phases. An im- 
portant point is that the amplitudes of the sine waves and 
their phases are constants of the motion, provided that the 
motion is linear. The example of Fig. 1, which has five sine 
waves, is a solution of equation (1) at time ¢ = 0. 


Cnoidal Waves and Nonlinear 
Fourier Analysis 


Unfortunately linear Fourier analysis has its limitations, 
primarily because shallow water waves are really not linear. 
The waves are long with respect to the depth and hence "feel" 
the presence of the bottom; this effect leads to the presence of 
nonlinear dynamics. One way that nonlinear effects manifest 
themselves ts by their influence on the shape of the waves. 
Generally speaking, shallow water waves have crests which 


are higher and narrower that those for sine waves; by the same 
token the troughs are less deep and broader relative to sine 
waves. For waves which are sufficiently high and steep break- 
ing can occur and consequently the free surface, momentarily, 
may not even be well defined; this is of course the situation in 
the surf zone. A complete understanding of nonlinear effects 
in surface water waves has proven to be an extremely difficult 
mathematical problem; it is enough to recall that the Euler 
equations for surface wave motion, marvellously complex in 
their nonlinearity, have remained unsolved for over two cen- 
turies. 

Since the Euler equations are so difficult to solve one often 
turns to simpler linear or nonlinear wave equations in order to 
more easily understand certain physical effects. In shallow 
water one of the most famous nonlinear equations is that due 
to Korteweg and deVries! ?; 


N+ CoN, + onn, + BN, ,,=0 (4) 


The constant coefficients here are the same as in equation 
(1) where a = 3 cp/2A. Note that equation (4) is the same as 


(1) except for the addition of the nonlinear term onn,. Of 


course the solution of (1) for periodic boundary conditions 
is trivial using the linear Fourier transform. The general IST 
solution to (4) for periodic boundary conditions required an 
additional 80 years of mathematical progress *''. A simple 
periodic solution of (4) is the cnoidal wave, which is the 
typical, classical nonlinear wave form often used in shallow 
water oceanography and offshore engineering “”: 


ata pacer = cos ( nkx) = 


n=) 


Qngen®{( K (m)/m)[ k(x- Ct) 3m} 


The modulus m of each elliptic function cn, the nonlinear 
phase speed C and the nome g are well-known formulas 
written in terms of the amplitude n, ” '; for brevity we do 
not give them here. The series representation for the cnoidal 
wave, given in (5), is the Stokes series solution to the KdV 
equation (4). When the modulus m — 0 the cnoidal wave 
reduces to a sine wave; when m —> 1 the cnoidal wave ap- 
proaches a solitary wave or soliton; intermediate values of 
the modulus correspond to the Stokes wave (see Fig. 2 for 
examples). Oceanographic applications of solitary waves 
were elegantly discussed in one of the early fundamental 
papers by Munk 2° 
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Figure 3. 


An example of nonlinear Fourier analysis using the inverse 
scattering transform. There are five cnoidal wave components 
in the spectrum. In (a) are the amplitudes and the moduli of 
the spectral components. In (b) are the cnoidal waves and 
their nonlinear interactions. Linear superposition of the cnoidal 
waves plus interactions yields the synthesized wave train at 
the bottom of the panel. 
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A key question at this point is the following: given that 
Fourier analysis for linear sine waves exists, does there also 
exist a nonlinear Fourier analysis for cnoidal waves? The 
major purpose of this paper is to discuss the fact that not only 
does nonlinear Fourier analysis exist, but it also has been 
formulated in a physical and mathematical form simple 
enough that practical oceanographic and engineering applica- 
tions of the approach are now under way. To show how this 
formulation arises I give the general solution to the KdV equation 
(4) in terms of the so-called 6-function representation 26 
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a2 
An (xt) = 2 5 In Oy (114.N2,---Mn)» 


where the 6-function is given by: 


Oy (1.N2---M) = 


> op 


N N N 
iM +> DL LM Binn Mn 
m=\in=1 


M, ...My=— n=1 


Here N is the number of degrees of freedom (or number of 
cnoidal waves) in a particular solution to the KdV equation. 
The 6-function phases are exactly the same as in linear 
Fourier analysis (see equation (3)): 


N, =k,xX— Ot + D, 


where the period matrix B = {B,,,,}, the wave numbers &,,, the 
frequencies @, and the phases @, depend upon algebraic 
geometric loop integrals whose determination is discussed 
elsewhere ’. The period matrix B is constant and provides the 
cnoidal wave amplitudes and their nonlinear pair-wise inter- 
actions. 

Note that the KdV equation (4) is a nonlinear generaliza- 
tion of the linear wave equation (1). Furthermore the cnoidal 
wave (5) is a nonlinear generalization of the sine wave (2). 
Finally, the 6-function formulation (6) is a nonlinear generali- 
zation of the linear Fourier series (3). 

It is straightforward to address the following theorem ’: 


Nonlinear Fourier Analysis Theorem: The 6-function solu- 
tion (6) to the KdV equation (4) can be written in the 


following form (u(x,f) = An(x,t), where A = o/68 = 3/2h%): 


a 
u (xt) =25 In @y (71) = lon (1) + Ming (TPB) 


Linear P Nonlinear 


cnoidal waves = among the cnoidal waves 





The basic result can be stated in words: Shallow water wave 
trains can be represented by a linear superposition of cnoi- 
dal waves plus a term which includes the mutual nonlinear 
interactions among the cnoidal waves. When the wave am- 
plitudes are small with respect to the depth, the cnoidal 
waves reduce to sine waves and the interactions tend to zero; 
consequently, linear Fourier analysis is recovered! It is in this 
sense that (8) is a nonlinear generalization of linear Fourier 
analysis. The cnoidal wave amplitudes are determined by the 





diagonal elements of the interaction matrix B while the 
interactions are determined by the off-diagonal terms. 

It is interesting to contrast the results of this Theorem with 
linear Fourier analysis. In this latter procedure the sine waves 
may be viewed as a set of linear basis functions. In the 
nonlinear problem the basis functions are no longer sine waves 
but are instead a nonlinear generalization, the cnoidal wave. 
According to the Theorem the nonlinear theory differs in two 
ways from ordinary linear Fourier analysis: it generalizes the 
basis functions from sine waves to cnoidal waves and includes 
the effect of the nonlinear interactions. This is a marvellously 
simple result. Fluid dynamicists, oceanographers and engi- 
neers are all familiar with the cnoidal wave. Now they can also 
do Fourier analysis with them thanks to inverse scattering 
theory and recent related advances ’. 

The remainder of this paper is devoted to the use of the 
above Theorem in both the construction of nonlinear wave 
trains and in the analysis of shallow-water ocean wave data. In 
the next Section I discuss a simple five-cnoidal wave solution 
in some detail; a twelve-cnoidal wave solution is also ad- 
dressed. Subsequently, I give the IST analysis of measured 
shallow water waves taken in the Adriatic Sea about 20 km 
from Venice, Italy. 





Figure 4. 


An example of two small-amplitude cnoidal waves from Fig. 
3 and their interactions. 
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Figure 5. 


An example of the interaction between the soliton and a sine 
wave of Fig.3. 
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Examples of the Synthesis of 
Nonlinear Wave Trains using 
Nonlinear Fourier Analysis 


The nonlinear Fourier synthesis of a wave train with 5 
cnoidal waves is given in Fig. 3. The cnoidal waves them- 
selves are shown vertically arranged in the upper half of the 
Figure. The uppermost component has a modulus which is 
near 1 and is therefore a solitary wave; the remaining com- 
ponents have intermediate values of the modulus and hence 
are Stokes waves. Below the cnoidal waves in Fig. 3 are 
shown the nonlinear interactions; the sum of the cnoidal 
waves plus the nonlinear interactions provides the synthe- 
sized wave train itself, as shown in the lower part of the 
figure. In what sense do the nonlinear interactions partici- 
pate in the dynamics of the wave train? It can be shown that 
the interactions have a very simple role: they phase shift the 
cnoidal waves in a collective fashion in order to construct a 
wave train which is an exact solution to the KdV equation. 
Thus the interactions shift each cnoidal wave to the right or 
the left as a function of space and time in order to include the 
nonlinear coupling among the cnoidal wave spectral compo- 
nents. 

There is a simple way to see how the nonlinear phase 
shifting occurs. Since the interactions are all pair-wise we can 
look at each pair of cnoidal waves and their nonlinear interac- 
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tions (see also Boyd 7’ and cited references). In Fig. 4 I give 
the two smallest components of Fig. 3; these waves are so 
small that they differ very little from ordinary sine waves. 
However, they do differ slightly from pure sine waves, since 
the nonlinear interactions, while small, are still significant. The 
subtle interplay between these two "small-amplitude Stokes 
waves" is quite nice and interesting. Summation of these waves 
plus interactions constructs their contribution to the 5 degree 
of freedom wave train of Fig. 3. 

Continuing this line of argument in Fig. 5, I take the first 
component of Fig. 3 (a soliton) and the last component (prac- 
tically speaking it is a sine wave) and compute the nonlinear 
interactions. This constitutes the interaction of a soliton and a 
sine wave. This has been accomplished by extracting the 
appropriate components from the B matrix for the example of 
Fig. 3 and then reapplying equations (6). Note the symmetric 
form of the nonlinear interactions about the peak of the soliton. 
Superposition of the soliton and the sine wave plus the non- 
linear interactions yield the synthesized wave train in the 
bottom of the figure. 

Fig. 6 shows the interaction of the first two components 
of the 5 cnoidal wave example of Fig. 3. Here the components 
are a soliton and a Stokes wave. The interactions are easily 
computed and the synthesized wave form is shown at the 
bottom of the figure. Finally a Stokes wave and a sine wave 
are studied in Fig. 7. The interactions are computed as before 
and the resultant synthesized wave train is also given. From 
these examples we see-that the interactions, as computed in 
pair-wise fashion, are typically stronger between larger cnoi- 
dal wave components. The more nonlinear the two cnoidal 
waves, the greater the pair-wise nonlinear interactions. 

Next I give the synthesis of a cnoidal wave train with 
12 components. The spectrum is shown in Fig. 8(a) where 
the cnoidal wave amplitudes (solid line) and the spectral 
moduli (dotted line) are graphed as a function of wave 
number. This graph differs qualitatively in only one way 
from that of the linear Fourier transform, e.g. the moduli 
indicate how nonlinear the individual spectral components 
are. In this example the first two components are quite 
nonlinear because their moduli are > 0.97; these compo- 
nents behave like solitons. The remaining components are 
substantially less nonlinear as the moduli fall below = 0.4; 
hence these components may be viewed as Stokes waves. 
The 12 cnoidal waves are shown in Fig. 8(b). The nonlinear 
interactions have also been computed; note that the largest 
pulses (effectively largest phase shifts) in the interaction 
wave train occur for the two peaks in the second cnoidal 
wave near the top of the figure. The final wave train is 
synthesized at the bottom of Fig. 8(b). 
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An example of the interaction between the soliton and the 
Stokes wave of Fig. 3. 
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Figure 7. 


An example of the interaction between a sine wave and the 
Stokes wave of Fig. 3. 
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Figure 8. 


An example of the nonlinear Fourier synthesis of a wave train 
with 12 cnoidal waves. In (a) are the spectrum and the moduli 
of the components as a ‘function of wave number. In (b) are 
shown the cnoidal waves, the nonlinear interactions and the 
synthesized wave train. 
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Figure 9. 


A measured surface wave time series from the Adriatic Sea. 
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Analysis of Shallow Water 
Ocean Waves with the Non- 
linear Fourier Technique 


In the last Section I discussed a number of simple exam- 
ples of the synthesis of nonlinear wave trains using cnoidal 
waves. However, a major goal of this research effort is to also 
analyze real ocean wave data using the IST technique. In Fig. 
9 I show a measured wave train from the Adriatic Sea in water 
of depth 16.5 m '*:!?. The significant wave height is 1.5 m; the 
largest waves are about 2.5 m in height. The average zero 
crossing period is 10.1 sec. The inverse scattering transform 
spectrum is given in Fig. 10. The cnoidal wave spectral ampli- 
tudes (solid line) and the moduli (dotted line) are shown. Note 
the behavior of the modulus, which is large at low frequency 
and then generally decreases with increasing frequency. The 
moduli also peak near the peak of the spectrum, rising to about 
0.2 at ~ 0.1 Hz. The interpretation of these results is that at 
low frequency one has several small amplitude solitons in the 
spectrum, while for the remainder of the spectrum the compo- 
nents can be viewed as Stokes waves of varying degrees of 
(relatively small) nonlinearity. Fig. 11 gives the first 50 cnoidal 
waves in the spectrum; truncating the number of cnoidal waves 
at fifty corresponds to a low pass filtering of the measured 
wave train in the frequency range 0 to 0.2 Hz. These cnoidal 
waves have been summed and the result is shown below the 
cnoidal waves in the figure. The nonlinear interactions have 
also been computed and these are shown below the summed 
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Figure 10. 


Nonlinear Fourier analysis of the Adriatic Sea time series of 
Fig. 9. Shown are the cnoidal wave spectrum (solid line) and 
moduli (dotted line) as a function of frequency. 
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cnoidal waves. Finally the low-pass-filtered time series has 
been reconstructed at the bottom of the figure. 

What have we learned from this exercise? It is clear that 
the nonlinearity in the measured wave train is not very great. 
This can be seen in two ways. First, in Fig. 10 one sees that 
the moduli are small, except for the first two components. A 
mild enhancement in the modulus occurs near the peak of the 
spectrum, a result which is to be expected for these larger 
waves. Second, the train of nonlinear interactions in Fig. 11 is 
rather small compared to the reconstructed wave train itself. 
However, even though the nonlinear effects are not very large, 
they are nevertheless substantial and provide some insight 
about the behavior of the measured wave train: Note that the 
sum of the cnoidal waves has larger amplitudes than does the 
wave train itself. This occurs because the interaction term is 
out of phase with the measured time series. This should not be 
surprising; we are all familiar with the two soliton interaction 
in which the solitons reduce their amplitudes during a colli- 
sion. As just seen this effect is also observed in the analysis of 
this Adriatic Sea wave train. An interesting observation is that 
because the interactions are relatively small in this case, one 
might eventually be able to construct a perturbative theory for 
the nonlinear interactions. 

What should happen as the waves propagate into shal- 
lower water? As the depth decreases the Ursell number 
(“Ve ) must rapidly increase. This of course increases the 
influence of the kottom on the nonlinear dynamics of the wave 
train. Consequently, on the average the moduli of the cnoidal 
waves must increase and consequently the nonlinear interac- 
tions will also increase substantially. Eventually at some depth 
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all the moduli will have reached 1 and all the cnoidal waves 
will have become solitons. Thus the term "nonlinear shoaling" 
takes on new meaning in the context of the inverse scattering 
transform. The problem of nonlinear shoaling has been pre- 
viously addressed, not for cnoidal waves, but for the hyperel- 
liptic function representation in references '*:!9, 

It goes without saying that IST provides new opportuni- 
ties for enhancing our knowledge and understanding of non- 
linear shallow water waves. Future efforts will address a whole 
host of interesting nonlinear effects, including resonances and 
chaos. 


Other Applications of Inverse 
Scattering Theory 





Figure 11. 


Nonlinear Fourier analysis of the Adriatic Sea time series of 
Fig. 9. Shown are the first fifty cnoidal waves in the spectrum; 
this is equivalent to a low pass filter of the measured wave data 
from 0 to 0.2 Hz. Also shown is the wave train corresponding 
to the sum of the cnoidal waves, the nonlinear interactions and 
the reconstructed low-pass-filtered input time series. 
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It is worth briefly mentioning that other applications of 
the inverse scattering technique are readily addressed in the 
shallow water coastal zone. These include the nonlinear dy- 
namics of internal waves, acoustic wave propagation in the 
shallow water sound channel and the nonlinear dynamics of 
surface and subsurface floating vessels. This discussion fo- 
cuses on present and future directions of research using IST. 

The shallow water internal wave field is often describable 
by the dynamics of the Korteweg-deVries equation 7*. The 
discovery of large internal solitons by Exxon drilling opera- 
tions in the Andaman Sea, offshore Thailand, has graphically 
illustrated the importance of internal waves on operational and 
design conditions for offshore operations (Osborne, Burch and 
Scarlet 7%). In this regard a soliton interpretation of the An- 
daman Sea data, using fundamental results of inverse scatter- 
ing theory, was made 15 years ago by Osborne and Burch **. 
Subsequently large internal solitons were also found in the 
Sulu Sea *° and other continental shelf areas. It thus seems 
timely that the inverse scattering transform be used to analyzed 
data of this type in order to enhance our understanding of the 
solitonic component of the internal wave field. 

Acoustic wave propagation on the continental shelf, in 
regions where KdV dynamics might be applicable, is also an 
important area of research. It is well-known that internal wave 
temperature variations in shallow water are often found to obey 
the Korteweg-deVries equation; because sound speed vari- 
ations closely follow the temperature variations it is clear that 
both temperatures and sound speed variations have inverse 
scattering transform modes. Therefore nonlinear Fourier 
analysis using the inverse scattering transform can provide 
additional perspective about sound propagation in these re- 
gions. 

Finally let me discuss floating body motions. How should 
the motion of a floating structure be influenced by nonlinear 
wave motion governed by IST? Can one compute the fre- 
quency response function of the vessel under the presence of 
wave forces due to a nonlinear wave train? The answers to 
these questions can be addressed in terms of the cnoidal wave 
modes and their mutual interactions. Let me give a simple 
illustrative example. Assume that the equation of motion of a 
floating object in a sea state governed by the KdV equation is 
given by 


mx+cx+kx=F(t) 


where F(t) is representable in terms of cnoidal waves plus 
interactions. Then the frequency response function of the 
floating body is easily determined in terms of the parameters 
m, c, k and the scattering transform modes. Therefore float- 
ing body motions, as governed by the above equation, may 
be described in terms of the nonlinear modes of the inverse 
scatting transform. Inverse scattering transform modes pro- 
vide a new and unique way to study important problems in 


ocean engineering. Many more vessel-motion problems are 
discussed in detail elsewhere *'. 


Summary and Discussion 


Given the rather simple physical interpretation of non- 
linear shallow water wave motion presented herein, why has 
inverse scattering taken so long to develop into data analysis 
tools? This is primarily because the development of the ap- 
proach has required a number of important new steps relating 
to the practical implementation of nonlinear time series analy- 
sis procedures. These steps include: 


(1) Closure of all aspects of the mathematical physics. 

(2) Development of the requisite numerical methods. 

(3) New and insightful experiments. 

(4) Development of fast time series analysis algorithms. 

(5) Development of physical insight with regard to the 
predictions of IST. 

(6) Development of procedures for the implementation of 
the algorithms, particularly with regard to nonlinear filtering 
and data analysis applications. 

Many of these topics and others are addressed in some of the 
cited references. 

Further directions in the future application of inverse 
scattering theory include: 

(1) The continued development of faster numerical meth- 
ods. 

(2) Study of resonances in nonlinear shallow water waves. 

(3) Study of the onset of chaotic motions. 

(4) Extension of the method to two dimensional wave 
motion and development of the concept of nonlinear direc- 
tional spectra. 

(5) The analysis of large quantities of data, including 
shallow water surface and internal waves, acoustic waves, 
floating vessel motions, etc. 

(6) Extension of the inverse scattering technique to higher 
order for the study of larger and shorter waves. 
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Effect of Wind and Shear on 


Surface Waves 


Maria Caponi, TRW Space and Technology Division, Redondo Beach, CA 
Philip Saffman, California Institute of Technology, Pasadena, CA 


Introduction 


“However, his (Tevake, a Polynesian navigator) most in- 
teresting concept to my mind is that the shape of waves in the 
open sea can sometimes indicate the presence and direction of 
a current.” 

Some of the primary sensors for remote investigation of 
the ocean properties are high frequency radars that respond 
only to wind sensitive short waves, in the wavelength range of 
1 to 10 cm. Hence, the interpretation of radar images of the sea 
surface requires the understanding of the generation of short 
waves by wind and their interaction with longer waves and 
currents. In most oceanic situations, the wind blows over 
pre-existing waves that are long enough to withstand dissipa- 
tion. The presence of these waves can significantly affect the 
interaction properties between the wind and the shorter wind- 
generated waves. Further, the distribution of the short waves 
on the phase of the long waves has a strong impact on the 
interpretation of high frequency radar imaging. 

The study of the effect of wind and shear on water waves 
can be divided into three areas. First, the analysis of the shape 
of waves of permanent form. Second, the investigation of the 
stability of a uniform flow with wind and shear to infinitesimal 
disturbances. Third, the stability of finite amplitude distur- 
bances as well as the stability of the finite amplitude waves of 
permanent form themselves. Of particular interest here is the 
effect of the presence of long waves of permanent form on the 
shape and stability of short waves. 


The classical theory of non-linear water waves can be said 
to have started with the work of Sir Gabriel Stokes in 1847, 
who calculated the first few terms in the Fourier series expan- 
sion of the wave surface for irrotational (i.e. vorticity free) 
gravity waves of permanent form, and conjectured that there 
would be a wave of greatest height which would be cusped. In 
the subsequent nearly 150 years, the field has attracted the 
attention of leading mathematicians and physicists, including 
e.g. Einstein”. The study of water wave phenomena has stimu- 
lated the discovery of novel methods for the exact solution of 
non-linear partial differential equations (e.g. soliton and in- 
verse scattering theory). The recurrence phenomenon of non- 
linear dynamical systems, known as the Fermi-Pasta-Ulam 
phenomenon has been found in water wave experiments and 
theory (e.g. Yuen and Lake’). Theoretical models of non-linear 
water wave interactions studied by Caponi, Saffman and Yuen* 
display a rich chaotic behavior. 

A spectacular experiment by Su et al. ° exhibited bifurca- 
tion of two-dimensional finite amplitude water waves into 
three dimensional waves of permanent form (figure 1). These 
waves, by a remarkable coincidence, were being calculated 
independently at the same time by Meiron, Saffman and Yuen® 
(figure 2). But note that the comparison between the experi- 
mental observations and the theory displays an important 
difference which is not yet adequately resolved; the theoretical 
waves have fore and aft symmetry, while this is not.so for the 
physical waves. It is an open question whether the symmetry , 
breaking is due to dissipation or nonlinearity. 
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With the exception of studies directed specifically at the 
generation of waves by wind, the classical theory generally 
ignores the dynamics of the air (because the air density is 
O(10°) smaller than that of the water) and neglects the possi- 
bility of shear in the water because the sources of vorticity are 
weak and the flow is irrotational to a good approximation, at 
least until wave breaking occurs. However, it is becoming 
recognized that there may be a number of important phenom- 
ena where the dynamical importance of the air and water shear 
plays a significant role. 

“A number of phenomena are still quite unexplained, and 
are very relevant to engineering problems. For example, it has 
been known by many generations of sailors that waves pro- 
duced by wind when it opposes a quite small current are 
greater than those when the same wind is in the same direction 
as the current. With ‘wind against tide’ the sea appears 
rougher, small ships toss about far more and there is more 
broken water seen. The effect is much more marked than can 
be accounted for by the slight increases of relative velocity. 
Sailors are taciturn folk, and so there are few references to this 
in the literature, and none at all in scientific work.’ Francis’. 

In the case of wind generated waves co-existing with a 
swell, an interesting phenomenon was experimentally demon- 
strated by Mitsuyasu®. The growth of the wind waves can be 
intensified or attenuated depending whether the swell is propa- 
gating against or in the direction of the wind. For the latter case 
and sufficient swell steepness (nonlinearity) the short wind 
generated waves and roughness are swept clean, pushed to the 
forward face of the steep long wave; see Mitsuyasu*, Hsu’. 

The principal theories that explain the transfer of energy 
from wind to waves have neglected until recently the effect of 
water shear and long waves. Part of the problem has been the 
need for analytical theories to use simplifying approximations 
that limit their range of application. For numerical solutions 
the difficulty is related to the extensive computations required 
for an accurate solution of the eigenvalue problems (i.e. wave 
speed calculations), even for the linear case if a continuous, 
realistic wind profile is assumed. This limits the possibility of 
obtaining a complete picture of the dominant effects. 

Recently, Caponi et al '° have analyzed the stability of the 
air-water interface using idealized wind and water shear profiles 
that reduce the linear eigenvalue problem to the solution of a 
quartic equation. The simplification also permits investigation of 
the effect of finite but smali steepness long waves on the overall 
stability of the short wind generated waves. The simplicity of the 
model allows a detailed parametric investigation that can be used 
to focus ideas and compare with detailed numerical investigations 
(e.g. Morland and Saffman"', Valenzuela '). 

The present article describes some of the current research 
and results associated with the effect of wind, shear and long 
waves on waves of permanent form and short wave instabili- 
ties. The description is concerned with inhomogeneity of 
velocity; effects of density stratification are not considered, 
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except for the separation into air and water both assumed of 
constant density. This is not because stratification is unimpor- 
tant; the contrary is true; but rather because we wish to focus 
here on effects of shear. 


Waves of Permanent Form in 
the Presence of Shear 


There have been many calculations of the shape of finite 
amplitude, two-dimensional waves of permanent form in the 
absence of shear and air motion (e.g. Cokelet'*). When the 
restoring force is gravity and surface tension is negligible, they 
are called Stokes waves. Calculations performed during the 
last 25 years have demonstrated that there is a continuous 
family of Stokes waves of given wavelength, from infinitesi- 
mal waves up to the wave of greatest height h/A = 0.14107, 
whose crests are cusped with interior angle 120°. h is the wave 
height, defined as the vertical distance between crest and 
trough, and A is the wavelength. 

The work of Chen and Saffman'* showed that the family 
of Stokes waves bifurcates subharmonically into waves of 
longer wavelength whose crests and troughs are not all of the 
same height and that the waves are not unique for h/A. > 0.1289. 
This is contrary to the long held belief that Stokes waves are 
unique (i.e. there is only one wave of given height and wave- 
length apart from a horizontal translation). It was also thought 
that the Stokes waves had fore and aft symmetry so that from 
the shape of a steady wave it could not be determined which 
way it is propagating. Zufiria!>, however, showed that symme- 
try breaking could occur spontaneously on the bifurcated 
branches so that waves without fore and aft symmetry exist. 





Figure 1, 


Bifurcation of two dimensional water waves into three di- 
mensional waves of permanent form. Experimental result. 
(From Su et al’). 











The amount of symmetry breaking is, however, quite small, 
being about 5% for the steepest wave that can be calculated 
before the accuracy of the solutions becomes doubtful. 

Less work has been done for the case when there is shear 
or vorticity in the water. For the case of constant vorticity and 
finite constant depth, for which the undisturbed velocity pro- 
file is a linear shear, numerical and analytical methods for 
small steepness waves were developed by several authors. The 
results from these studies did not show significant differences 
from those obtained without shear. 

The first work for the case of infinite depth and large 
steepness waves was carried out by Simmen and Saffman’®. 
They used a boundary integral method for a linear shear of 
infinite extent. This method expresses the solution in terms of 
an integral of the values of the fluid velocity at the boundary 
and leads to an integral equation that must be solved numeri- 
cally. Milinazzo and Saffman'’ used a Fourier series repre- 
sentation of the surface displacement for the case of a linear 
shear of finite thickness above infinitely deep water, the object 
being to model a wind drift layer (i.e., a shallow current layer 
in the ocean ostensibly produced by the drag of the wind). The 
interesting result was found that shear steepens waves propa- 
gating in the direction opposite to the drift current and en- 
hances the tendency for the wave to break. On the other hand 
the wave is flatter when drift current and propagation speed 
are in the same direction and the tendency to break is retarded. 
This effect can be quite strong. The limiting wave height can 
be reduced to one third of the no shear case when the wind drift 
velocity is around 50% of the wave speed, and increased by a 
factor of three when the wave and drift directions are opposite. 
The results also confirm a heuristic one-dimensional model of 
breaking by Banner and Phillips'*. A formulation using a 
boundary integral method was also carried out for finite depths 
and large steepness by Teles da Silva and Peregrine'’. The 
results obtained from this study are similar to those obtained 
from the infinite depth case. 

The families of waves of permanent form with shear also 
exhibit subharmonic bifurcation into longer waves with more 
structure. Indeed, because the shear makes the undisturbed flow 
non-symmetric, it could be anticipated that the breaking of fore 
and aft symmetry occurs in the main branch family (i.e. waves 
that are analytic continuations with respect to shear of the Stokes 
waves). Surprisingly, this is not the case and the main branch 
family has fore and aft symmetry for all heights up to the 
maximum. Symmetry breaking does exist along the subharmonic 
bifurcation branches and is enhanced by shear (Baumstcin™. 
Figure 3 shows a main branch wave and a subharmonic symmetry 
breaking member of the bifurcated branch. 

To date, the calculations for the shapes of waves of perma- 
nent form are restricted to two-dimensional disturbances of 
piecewise linear, undisturbed velocity profiles (sometimes called 
stick profiles) for mathematical reasons. In this case, vorticity is 
a conserved quantity, which does not need to be solved for, and 


the disturbance is irrotational. The mathematical techniques of 
potential theory can then be applied and the problem formu- 
lated using complex variable theory. In particular, boundary 
integral methods can be employed, that are free of the conver- 
gence difficulties that can arise with the use of truncated 
Fourier expansions. The choice of numerical method plays an 
important role in this type of calculation. However, it is not 
essential to restrict attention to stick profiles and current work 
by us is studying the form of permanent waves on continuous 
profiles of velocity or vorticity. The extension to three-dimen- 
sional waves is a very difficult problem which appears to be 
beyond present analytical and numerical capabilities. 


Effect of wind and water shear 
on wave instabilities 


There are two aspects of the way in which wind and shear, 
or vorticity in the air and water, affect the instabilities of 
surface waves at an air-water interface. 

First, there is the question of the instability of existing 
waves, generated mechanically for example, or by earlier 
winds. In this situation, the dynamics of the air is unimportant 
and the hydrodynamics of the water motion is the dominant 
factor. The second aspect is concerned with wave generation 
by wind or current. In particular, we can consider the genera- 
tion of waves as a problem of the linear hydrodynamic insta- 
bility of a laminar flow. 


Instabilities of existing waves 


The instadilities of wave trains can be categorized, for 
convenience, into three types, as follows. 





Figure 2. 


Bifurcation of two dimensional water waves into three dimen- 
sional waves of permanent form. Theoretical result. (From 
Meiron, Saffman and Yuen®) 
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(i) Resonant interactions, Phillips”': This describes the 
transfer of energy between wave trains of vector wavenumbers 
k, and frequencies @, when )'k, = 0 and }\@, = 0. 

(ii) Two- and three-dimensional long wave envelope 
modulations or side band instabilities: These go generally 
under the name of the Benjamin and Feir” instability. Insta- 
bilities of this kind were also described by Lighthill” (who 
used Whitham’s variational principle), Zakharov’* and oth- 
ers. 

(iii) Two- and three-dimensional, sub- and super-har- 
monic instabilities of finite amplitude waves of permanent 
form: The two-dimensional case was studied by Longuet-Hig- 
gins”. Saffman”’ and MacKay and Saffman”* applied a 
Hamiltonian formulation to explain qualitatively some of the 
results such as the change of stability when the wave energy 
as a function of wave height is stationary or modes collide. 
This type of formulation makes it possible to follow a change 
in the eigenvalues (i.e. growth rate of the disturbance) as a 
control parameter such as the wave height is varied. The 
three-dimensional case was studied by McLean, Ma, Martin, 





Figure 3. 


Waves of permanent form with shear. a) Main branch wave. b) 
Subharmonic bifurcation branch illustrates symmetry break- 
ing due to shear. 
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Saffman and Yuen”’, see also McLean*®; we shall refer to this 
as the MMMSY instability. 

Although the phenomena in (ii) and (iii) can be regarded 
as particular cases of (i), it is useful to categorize them sepa- 
rately since for those cases the small disturbances grow at slow 
but exponential rates rather than the linear rate usually associ- 
ated with the resonances of (i). 

The resonances and envelope modulations are basically 
slowly varying, weakly nonlinear phenomena, and the modi- 
fications due to shear will arise primarily through changes in 
the linear dispersion relations (i.e. the relations between wave 
frequency and wavelength) and the speed of weakly nonlinear, 
uniform wave trains. For example, Li, Hui and Donelan*' have 
calculated the effect of constant water shear on the side-band 
instability for the case of infinite depth, and found that small 
shear tends to enhance instability whereas large shear tends to 
suppress it. 

The effect of shear on the MMMSY wave instabilities 
is a more difficult problem as potential theory cannot be 
used to relate perturbations of the flow to perturbations of 
the surface. Also, we cannot use the mathematical simplifi- 
cations due to existence of a scalar stream-function as it 
does not exist for three-dimensional long wave distur- 
bances. However, Benney and Chow™ have used perturba- 
tion methods for weakly nonlinear waves and found that 
weak shear generates relatively strong three-dimensional 
long wave instabilities. The properties of the instability are 
significantly different from those of the Benjamin and Feir 
type. For finite amplitude waves and zero surface tension, 
Okamura and Oikawa* have calculated the instability of 
two dimensional waves of permanent form to three-dimen- 
sional disturbances. They recover the finite depth MMMSY 
results (McLean™) as the shear tends to zero. They also see 
the Benney-Chow long oblique waves, which are absent 
when there is no shear. 


Wave generation by wind or current 


Three basic theories have been used until recently to 
describe the generation of waves by wind or water shear. 

(i) Kelvin-Helmholtz instability: This is the classic insta- 
bility of a vortex sheet, modified by a density difference 
between the fluids on the two sides of the sheet. 

(ii) Miles*® instability: This is the transfer of energy from 
a shear flow in the air over water initially at rest into waves 
whose phase speed lies between the minimum and maximum 
of the wind velocity. It has been investigated by Morland and 
Saffman'" as an inviscid hydrodynamic instability of a smooth 
profile. This formulation leads to an eigenvalue problem for 
the wave speed and growth rate involving a solution of the 
Rayleigh equation. Moder computers and methods which 
locate the eigenvalues in the complex plane make the problem 
easy to solve numerically and render unnecessary the hard 
analysis of the critical layer, where the fluid velocity equals 





the wavespeed. This is also the case when an exponential 
profile is chosen so that the Rayleigh equation has exact 
solutions in closed form in terms of hypergeometric functions. 

(iii) Drift layer instability: This is the hydrodynamic in- 
stability of a current with shear in the water, which could be 
produced, for example, by the drag of the air over the water or 
mechanically by running a stream over the surface, as for 
example in an estuary. It was investigated by Esch*® for gravity 
waves on piecewise linear profiles, by Stern and Adam*” for 
capillary-gravity waves on piecewise linear profiles and by 
Morland, Saffman and Yuen*® for capillary-gravity waves on 
smooth profiles. These treatments are all inviscid and the direct 
effect of viscosity on the wave evolution is neglected. Viscosity 
may, of course, play an important role in the generation of the 
shear profiles, but this is assumed to be on a much longer time 
scale than that of the instability. Viscosity is sometimes em- 
ployed to circumvent analytical difficulties associated with a 
critical layer; this introduces the Orr-Sommerfeld equation 
(e.g. Valenzuela’? and subsequent studies of coupled air-water 
shear flows). But a consistent purely inviscid treatment is 
possible and the critical layer need not be considered, as 
mentioned above. It is interesting that a necessary condition 
for the existence of the drift layer instability is that the drift 
layer surface velocity U,, say, is larger than c,, = (4gT)', the 
minimum phase speed of capillary-gravity waves (g is the 
acceleration due to gravity and T is the surface tension). It is 
also found that co-flowing waves are stable, and that the 
unstable waves tend to be moving very slowly. 

The mechanism for the generation of waves by drift layer 
instability is qualitatively attractive. It also provides a delay 
mechanism for the wave generation after a wind starts blowing, 
because it is found that the instability mechanism will not work 
until the layer thickness A reaches a critical value as shown by 
Morland, Saffman and Yuen** and Caponi et al**. But quanti- 
tative predictions are very poor. The necessary condition U/, > 
C,, IS not satisfied by the wind drift layer, and the predicted 
phase speeds are of order '4c,, which is much too fast to be 
consistent with observation. On the other hand, the theory may 
be applicable to tidal flow in the presence of wind. An onshore 
wind, for example, would tend to amplify counter flowing 
waves on an ebbing tide which would be unstable and produce 
choppiness, whereas unstable waves would not be supported 
with an onshore wind and a rising tide and the surface would 
remain relatively smooth. 

A further mechanism (Phillips*°), in which waves are 
generated by random pressure fluctuations due to a turbulent 
wind will not be considerer here. For this case, the main effect 
of shear is the small modification of the dispersion relation. 

Recently a model has been developed that encompasses 
the three main basic theories and because of its simplicity 
allows a detailed parametric investigation of the characteristics 
of wind generated waves, adding new insights of benefit to 
remote sensing. This model is described in the next section. 


A simple model of wind 
shear interaction 


The difficulty of the problem is increased when both wind 
and water shear are included in the analysis of the generation 
of water waves. A model by Caponi, Caponi, Saffman and 
Yuen!®, hereafter called the CCSY model, makes use of a 
piecewise linear profile to simplify the problem and has led to 
the prediction of interesting new phenomena in the behavior 
of wind generated water waves. The advantage of using such 
a profile is that the system dispersion relation becomes a 
quartic equation in the complex frequency and the parametric 
dependence of the solutions can be easily examined. Although 
models that make use of smooth profiles can result in more 
accurate growth rates they require at best (e.g. for the case of 
exponential profiles) the solution of dispersion relations con- 
sisting of transcendental equations involving hypergeometric 
functions. An overall understanding of the effects involved 
becomes more limited the more complex the dispersion equa- 
tion. This is due to the extensive computations necessary to 
obtain a detailed dependence of the wave characteristics (e.g. 
phase velocity, growth rate, wavelength) on the various pa- 
rameters of interest (e.g. shear layer thickness, wind velocity, 
etc.) 

The idealized undisturbed horizontal velocity profile used 
in the CCSY model is (see figure 4) 


0 if y<—D,,; 

U,, + Uwy/,, if -D,<y<0 

U,,+Uay/, if 0<y<D, 
[U,,+U, 


U(y) = 


ify>D, 


This represents a double shear layer in the air and in the 
water. An infinitesimal wave-like disturbance with frequency 
@ and wave number k = 2 1/A is imposed, and an inviscid fluid 
assumed. Linearization about a flat surface reduces the mo- 
mentum and continuity equations to a single equation for the 
y-dependence of the vertical velocity v(y). In each of the four 
regions, the solution takes the form v,(y) = A,e” + B, e®, 
where n = 1,2,3,4 according to the region. The boundary 
conditions: continuity of the vertical displacement and pres- 
sure across all boundaries, and vanishing velocity at great 
height and depth are used afterwards to obtain a dispersion 
relation between the frequency @ and wave number k. The 
relation can be expressed as a quartic equation for the ratio of 
the complex wave frequency (real frequency + i growth rate) 
to the gravity—capillary waves frequency as a function of 
known constants (gravitational acceleration, surface tension, 
density of air and water), the relative water to air velocity and 
dimensionless parameters kD,, KD,, and the ratio of the wind 
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speed (U,) to the phase velocity of the gravity—capillary waves 
(cg). The assumption of continuity of tangential stress and a 
fixed percentage for the relative water to air velocity, U, = 
U,JU,, is used to reduce the results toa single stability diagram 
in kD, vs U,/co space. 

The results are illustrated in figures 5, 6 and 7. These 
figures show contours of constant y /@,, the ratio of the growth 
rate to the intrinsic frequency of the gravity-capillary waves 
at intervals of 0.1. Figure 5 corresponds to a water—air ratio of 
0 04, consistent with the common estimate derived from ob- 
servations of drift velocities. Figures 6 and 7 illustrate the 
variations in the stability diagram due to either larger or 
negligible U,. In the figures the instability region (the shaded 
area) is defined by y/@, # 0. 

The overall stability is strongly dependent on water drift 
and wind velocity. For example, in the absence of water drift 
(figure 7), the instability region generally comprises two re- 
gions arising from two separate pairs of complex roots of the 
quartic equation. Region A is consistent with the critical layer 
analysis of Miles *°. It occurs for winds below 6 m/s, the phase 
velocity of the unstable waves is approximately the phase 
velocity of the gravity capillary waves, the growth rate in- 
creases with U,/c, and the instability region narrows and 
asymptotes to U,/c, when kD, >>I. Region B represents an 
extension of the classical Kelvin-Helmholtz instability 
mechanism to finite drift layers. It occurs for winds above 6 
m/s, the wave’s phase velocity in this region is much smaller 
than cy and the growth rate much larger than in Region A with 
a parameter dependence in agreement with Kelvin—Helmholtz 
instability. 

In the presence of water drift (see figures 5 and 6) an 
intermediate region develops for winds between approxi- 
mately 4 m/s and 6 m/s, where the growth rate decreases as U, 
increases. Of particular interest is the existence of a range of 





Figure 4. 


Idealized wind profile and geometry for CCSY model. 
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Figure 5. 


Unstable region for water to air velocity ratio of .04 (shaded 
area). Solid lines are contours of constant ratio of growth rate 
to intrinsic frequency of the gravity-capillary waves at intervals 
of . 1. 























stable wavenumbers for U, > 0.03 and the deduction that 
different wave components of a disturbance behave differently 
under a given wind in the presence of water shear. This 
prediction can be critically tested by experiment. If validated, 
this prediction may suggest a new design for radar scatterome- 
ters with simultaneous multifrequency radars which can more 
accurately deduce the wind speed by taking advantage of the 
difference in behavior and history of backscatter returns from 
different frequency channels. The existence of regions where 
short waves surrounding the wave of minimum speed do not 
grow also indicates a potential effect of the presence of long 
waves on the short waves (cf. next section). 

A study similar to the CCSY analysis, that makes use of 
n multilayers and includes nonlinear effects, has been carried 
out by Gertsenshtein, Romashova and Chernyavskii*'. The 
results presented in that paper do not agree with those of CCSY 
and unfortunately the article does not present sufficient detail 
to understand the source of disagreement. 

A study of the instability for continuous profiles has been 
carried out by Hill*?. Apart from a change of scale of the 
growth rates, the results are qualitatively similar to those of 
CCSY. In this case the source of difference in growth rate is 





the infinite curvature assumed in the piece-wise linear profiles 
at a change of slope. 


Non linear interactions 


An interesting and important question is the alteration of 
the growth rate of the linear instabilities caused by wind and 
water shear when a long wave of permanent form is present. 
In the absence of wind and for stable long waves this is the 
modification of the MMMSY problem by non-uniform water 
shear. For capillary gravity waves, it is possible that there 
might be enhanced instability because water shear causes extra 
corrugation of the interfaces (Milinazzo and Saffman!”). Pre- 
liminary calculations indicate that the swell can reduce the 
growth rate of the MMMSY instability, but the effect is, 
however, not large. 

For the case of wind but no water shear, an extension of 
the CCSY analysis has been carried out to study the effect of 
pre-existing waves of wavelength A, on the stability of shorter 
wind generated waves of wavelength A, = A,/2 (Caponi *°). 
The analysis assumes infinitesimal wave perturbations riding 
on longer wave disturbances of finite amplitude a. The dis- 
placements at each interface i (see figure 4 ) are then assumed 
of the form),;=,+N, with N,>n,, where 
Nj = 4 COs (k,x). 


The dynamic and kinematic boundary conditions are ex- 
panded about the long wave displacement and only terms of 
order k,a are kept. A system of equations is obtained that is 
linear in the perturbations and has coefficients of periodicity 
2n/k, . These equations can be solved by expanding the pertur- 
bations in Fourier series of the form:  ,, = 51) ,,¢ *e” with 
N n41,, << 11 ,,;- These assumptions reduce the problem of find- 
ing the system dispersion relation to that of finding the zeroes 
of the determinant of a matrix M, whose coefficients are 
functions of the wave frequency, growth, wavelength and the 
system parameters. Unfortunately, for a system consisting of 
4 interfaces and keeping terms of order n in the expansion, M 
is of order (9 n). The sparcity of the matrix determinant results 
in a dispersion equation of order 4n in the complex frequency, 
which defeats the purpose of a simple model for large n. 

For a system with only wind shear and n=2, M is of order 
12 x 12, which is analytically manageable although algebrai- 
cally cumbersome. The problem can be reduced to a dispersion 
equation of order 6 in the complex frequency that has an 
interesting analytical form. Basically it corresponds to the 
linear (either CCSY or gravity capillary) dispersion relation of 
the waves of wavelength k, coupled to the linear (CCSY) 
dispersion relation of the waves of wavelength 2k, = k, by a 
term proportional to the longer wave slope squared. The cou- 
pling term is a function of the system parameters and the wave 





Figure 6. 


Unstable region for water to air velocity ratio of .057 (shaded 
area). Solid lines are contours of constant ratio of growth rate to 
intrinsic frequency of the gravity-capillary waves at intervals of . 1. 


Figure 7. 


Unstable region for water to air velocity ratio of 0.0 (shaded area). 
Solid lines are contours of constant ratio of growth rate to intrinsic 
frequency of the gravity-capillary waves at intervals of . 1. 
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Figure 8. 


Unstable region for along wave of ka = .02 moving in the same 
direction as the wind, water to air velocity ratio of 0.0 (shaded 
area). Solid lines are contours of constant ratio of growth rate 
to intrinsic frequency of the gravity-capillary waves at intervals 
of .1. Dashed lines enclose unstable region in the absence of 
the long wave (see figure 7). 





Figure 9. 


Unstable region for a long wave of ka = .02 moving in the 
opposite direction to the wind, water to air velocity ratio of 0.0 
(shaded area). Solid lines are contours of constant ratio of 
growth rate to intrinsic frequency of the gravity-capillary waves 
at intervals of .1. Dashed lines enclose unstable region in the 
absence of the long wave (see figure 7). 
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characteristics, its sign determines if the presence of the long 
wave enhances or reduces the growth of the shorter wave. The 
stability diagrams of figures 8 and 9 illustrate this effect. In 
these figures, the dashed lines enclose the unstable region in 
the absence of long waves, see figure 7. 

Figure 8 corresponds to a long wave with a slope of 0.02 
moving in the same direction as the wind. The effect is stabi- 
lizing, the instability region becomes smaller and the value of 
the growth rate inside the region is reduced. The effect is 
stronger the larger the value of the longer wave slope, the 
unstable region becomes negligible as k,a — 0.04. Figure 9 
illustrates the case of the effect of a long wave moving opposite 
the direction of the wind. In this case the growth rate and the 
growth region are enhanced. These results are very important 
because they offer a potential explanation for the sweep up 
mechanism (Yuen “*). Further, they are in agreement with 
recent experimental measurements by Mitsuyasu® and Hsu’. 
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For n=2, the Caponi analysis is currently being further ex- 
tended to include water shear in order to analyze the effect of 
long waves on the predicted stable region. 

For the case of arbitrary n it is possible to formally analyze 
the system to obtain the general form of the dispersion relation 
(Caponi “°). It was found that the general dispersion relation 
of the system could be written as a coupling between the linear 
(CCSY) dispersion relations for each harmonic. The coupling 
term in this case can be written as a sum of n terms each one 
proportional to the slope of the longer wave to the j power with 
j running from 2 to n. This result illustrates the coupling of 
higher harmonics with the fundamental long wave as well as 
the dependence on the long wave slope. The analysis, however 
is only formal and can not be used to obtain a quantitative 
result. 

In order to investigate the stability of short waves with 
A, << A, it is necessary to use a different approach. Troitskaya 





utilized a multiple spatial scale expansion to study the case of 
wind but no water shear. Her results seem to indicate that the 
growth rates associated with the Miles instabilities can be 
significantly reduced in the presence of the longer waves, again 
in agreement with observation. Caponi and Saffman * are 
currently investigating an expansion that makes use of the 
assumption that the effect of the very long waves on the 
unstable short waves occurs on a time scale that is slow 
compared to the short wave inverse frequency. The approach 
includes the effect of nonlinear wave-wave interaction and 
modulational instability promising to encompass previous 
nonlinear treatments. 
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RADM Gaffney 
Is New Chief of 
Naval Research 


Rear Admiral Paul G. Gaffney, II, became the 19th Chief of Naval Research on July 12, 1996, during a change of command at 
the Washington Navy Yard. He relieved Rear Admiral Marc Pelaez. 


In addition to his new position, RADM Gaffney will continue as Commander of the Naval Meteorology and Oceanography 
Command for the next year. He has held this position since August 1994. He is dividing his time between Office of Naval 
Research (ONR) Headquarters in Arlington, Virginia, and the Naval Meteorology and Oceanography Command Headquarters 
in Stennis Space Center, Mississippi. 


As the Department of the Navy’s senior science and technology executive, RADM Gaffney is responsible for more than a $1 
billion science program ranging from basic research, conducted at universities to development of manufacturing technologies. 
In addition to headquarters, RADM Gaffney’s ONR command includes offices in London and Tokyo which provide liaison to 
the international scientific establishment. Also under his command is the Navy’s corporate laboratory, the Naval Research 
Laboratory (NRL), in Washington, D.C., with branches in Stennis Space Center, Mississippi, and Monterey, California. NRL is 
a leader in the development of science and technology relevant to all aspects of naval warfare, including control of the seas, 
littoral areas, naval airspace, the middle atmosphere, and outer space. 


Celebrating 50 years of service to the Department of the Navy and the nation, ONR in 1996 meets the unique requirements of the 
Navy and Marine Corps by supporting research at Navy laboratories and funding research through grants and contracts with 
universities, industry, and non-profit organizations. ONR provides technical advice to the Chief of Naval Operations and the 
Secretary of the Navy, works with industry to improve technology and manufacturing processes to reduce fleet costs, and 
encourages academic interest in naval-relevant science from the high school through post-doctoral levels. 


RADM Gaffney is no stranger to ONR. He was the Acting Director of the Arctic and Earth Sciences Division at ONR from 
November 1980 to August 1981. From May 1989 to February 1991, he served as the Assistant Chief of Naval Research and then 
as Commanding Officer of NRL. His other assignments have included: tours as Military Assistant to the Assistant Secretary of 
Defense for International Security Affairs, Commanding Officer of Oceanographic Unit Four survey in Indonesia and Com- 
manding Officer, Naval Oceanography Command Facility Jacksonville, Director of the Oceanography Resources Division on 
the staff of the Oceanographer of the Navy, and Executive Assistant to the Oceanographer of the Navy. 


RADM Gaffney is a 1968 graduate of the U.S. Naval Academy. He received a masters degree in ocean engineering from 
Catholic University of America, Washington, D.C., and a master’s degree in business administration from Jacksonville Univer- 
sity, Florida. He graduated with highest distinction from the Naval War College in Newport, Rhode Island in 1979. His awards 
include the Defense Superior Service Medal, Legion of Merit (three awards), Bronze Star with Combat “V,” Meritorious Service 
Medal and the J. William Middendorf Prize for Strategic Research at the Naval War College. 




















